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SPATIAL AND SPECTRAL WAVEFRONT ANALYSIS AND MEASUREMENT 

REFERENCE TO RELATED APPLICATIONS 
This application is based on provisional application U.S. Serial No. 
60/196,862, filed on April 12, 2000. 

FIELD OF THE INVENTION 
The present invention relates to wavefront analysis generally and to various 
applications of wavefront analysis. 

BACKGROUND OF THE INVENTION 
The following patents and publications are believed to represent the current 
state of the art: 
US PATENTS: 

5,969,855; 5,969,853; 5,936,253; 5,870,191; 5,814,815; 5,751,475; 5,619,372; 
5,600,440; 5,471,303; 5,446,540; 5,235,587; 4,407,569; 4,190,366; 
NON-US PATENTS: 

JP 9230247 (Abstract); JP 9179029 (Abstract); JP 8094936 (Abstract); JP 7261089 
(Abstract); JP 7225341 (Abstract); JP 6186504 (Abstract); 

OTHER PUBLICATIONS: 

Phillion D.W. "General methods for generating phase-shifting interferometry 
algorithms" - Applied Optics, Vol. 36, 8098 (1997). 

Pluta M. "Stray-light problem in phase contrast microscopy or toward highly sensitive 

phase contrast devices: a review"- Optical Engineering, Vol. 32, 3199 (1993). 

Noda T. and Kawata S. "Separation of phase and absorption images in phase-contrast 

microscopy"- Journal of the Optical Society of America A, Vol. 9., 924 (1992). 

Creath K. "Phase measurement interferometry techniques" - Progress in Optics XXVI, 

348 (1988). 

Greivenkamp J.E. "Generalized data reduction for heterodyne interferometry" - Optical 
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Engineering, Vol. 23, 350 (1984). 

Morgan C.J. "Least- squares estimation in phase-measurement interferometry" - Optics 
Letters, Vol. 7,368 (1982). 

Golden LJ. ''Zernike test. 1: Analytical aspects"- Applied Optics, Vol. 16, 205 (1977). 
Bruning J.H. et aL "Digital wavefront measuring interferometer for testing optical 
surfaces and lenses" - Applied Optics, VoL 13, 2693 (1974). 

SUMMARY OF THE INVENTION 

The present invention seeks to provide methodologies and systems for 
wavefront analysis as well as for surface mapping, phase change analysis, spectral 
analysis, object inspection, stored data retrieval, three-dimensional; imaging and other 
suitable applications utilizing wavefront analysis. 

There is thus provided in accordance with a preferred embodiment of the 
present invention a method of wavefront analysis. The method includes obtaining a 
plurality of differently phase changed transformed wavefronts corresponding to a 
wavefront being analyzed which has an amplitude and a phase, obtaining a plurality of 
intensity maps of the plurality of phase changed transformed wavefronts and employing 
the plurality of intensity maps to obtain an output indicating the amplitude and phase of 
the wavefront being analyzed. 

There is also provided in accordance with a preferred embodiment of the 
present an apparatus for wavefront analysis including a wavefront transformer operating 
to provide a plurality of differently phase changed transformed wavefronts 
corresponding to a wavefront being analyzed which has an amplitude and a phase, an 
intensity map generator operating to provide a plurality of intensity maps of the 
plurality of phase changed transformed wavefronts and an intensity map utilizer, 
employing the plurality of intensity maps for providing an output indicating the 
amplitude and phase of the wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of intensity maps are employed to provide an analytical output indicating the 
amplitude and phase. 

Still further in accordance with a preferred embodiment of the present 
invention the plurality of differently phase changed transformed wavefronts are 
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obtained by interference of the wavefront being analyzed along a coraraon optical path. 

Additionally in accordance with a preferred embodiment of the present 
invention the plurality of differently phase changed transformed wavefronts are realized 
in a manner substantially different from performing a delta- function phase change to the 
transformed wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of intensity maps are employed to obtain an output indicating the phase which 
is substantially free from halo and shading off distortions. 

Preferably, the plurality of differently phase changed transformed wavefronts 
include a plurality of wavefronts resulting from at least one of application of spatial 
phase changes to a transformed wavefront and transforming of a wavefront following 
application of spatial phase changes thereto. 

Additionally in accordance with a preferred embodiment of the present 
invention, the step of obtaining a plurality of differently phase changed transformed 
wavefronts includes applying a transform to the wavefront being analyzed thereby to 
obtain a transformed wavefront and applying a plurality of different phase changes to 
the transformed wavefront, thereby to obtain a plurality of differently phase changed 
transformed wavefronts. Preferably, the plurality of different phase changes includes 
spatial phase changes and the plurality of different spatial phase changes are effected by 
applying a time- varying spatial phase change to part of the transformed wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of different spatial phase changes are effected by applying a spatially uniform, 
time-varying spatial phase change to part of the transformed wavefront. Preferably, the 
transform applied to the wavefront being analyzed is a Fourier transform and wherein 
the step of obtaining a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts includes applying a Fourier transform to the plurality of 
differently phase changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the wavefront being analyzed is a Fourier transform and the 
plurality of different spatial phase changes includes at least three different phase 
changes. Preferably, the plurality of intensity maps includes at least three intensity maps 
and the step of employing the plurality of intensity maps to obtain an output indicating 
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the amplitude and phase of the wavefront being analyzed includes: expressing the 
wavefront being analyzed as a first complex function which has an amplitude and phase 
identical to the amplitude and phase of the wavefront being analyzed, expressing the 
plurality of intensity maps as a function of the first complex function and of a spatial 
funcfion governing the spatially uniform, time-varying spatial phase change, defining a 
second complex function, having an absolute value and a phase, as a convolution of the 
first complex function and of a Fourier transform of the spatial function governing the 
spatially uniform, time-varying spatial phase change. Expressing each of the plurality of 
intensity maps as a third function of: the amplitude of the wavefront being analyzed, the 
absolute value of the second complex function, a difference between the phase of the 
wavefront being analyzed and the phase of the second complex function and a known 
phase delay produced by one of the at least three different phase changes corresponding 
to one of the at least three intensity maps, solving the third function to obtain the 
amplitude of the wavefront being analyzed, the absolute value of the second complex 
fi^nction and the difference between the phase of the wavefront being analyzed and the 
phase of the second complex function, solving the second complex function to obtain 
the phase of the second complex function and obtaining the phase of the wavefront 
being analyzed by adding the phase of the second complex function to the difference 
between the phase of the wavefront being analyzed and the phase of the second complex 
function. 

Further in accordance with a preferred embodiment of the present invention the 
absolute value of the second complex function is obtained by approximating the 
absolute value to a polynomial of a given degree. 

Still further in accordance with a preferred embodiment of the present 
invention the second complex function is obtained by expressing the second complex 
function as an eigen-value problem where the complex function is an eigen-vector 
obtained by an iterative process. 

Preferably the second complex function is obtained by: approximating the 
Fourier transform of the spatial function governing the spatially uniform, time-varying 
spatial phase change to a polynomial and approximating the second complex function to 
a polynomial. 

Preferably, the wavefront being analyzed^ the absolute value of the second 
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complex function, and the difference between the phase of the second complex function 
and the phase of the wavefront being analyzed, are obtained by a least-square method, 
which has increased accuracy as the number of the plurality of intensity maps increases. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of different phase changes includes at least four different phase changes, the 
plurality of intensity maps includes at least four intensity maps and employing the 
plurality of intensity maps to obtain an output indicating the amplitude and phase of the 
wavefront being analyzed and includes: expressing each of the plurality of intensity 
maps as a third function of: the amplitude of the wavefront being analyzed, the absolute 
value of the second complex function, a difference between the phase of the wavefront 
being analyzed and the phase of the second complex function, a known phase delay 
produced by one of the at least four different phase changes in which each corresponds 
to one of the at least four intensity maps and at least one additional unknown relating to 
the wavefront analysis, where the number of the at least one additional unknown is no 
greater than the number by which the plurality intensity maps exceeds three and solving 
the third function to obtain the amplitude of the wavefront being analyzed, the absolute 
value of the second complex function, the difference between the phase of the 
wavefront being analyzed and the phase of the second complex function and the 
additional unknown. 

Further in accordance with a preferred embodiment of the present invention the 
phase changes are chosen as to maximize contrast in the intensity maps and to minimize 
effects of noise on the phase of the wavefront being analyzed. 

Preferably, expressing each of the plurality of intensity maps as a third function 
of: the amplitude of the wavefront being analyzed, the absolute value of the second 
complex function, a difference between the phase of the wavefront being analyzed and 
the phase of the second complex function and a known phase delay produced by one of 
the at least three different phase changes which corresponds to one of the at least three 
intensity maps and includes: defining fourth, fifth and sixth complex functions, none of 
which being a function of any of the plurality of intensity maps or of the time-varying 
spatial phase change, each of the fourth, fifth and sixth complex functions being a 
function of: the amplitude of the wavefront being analyzed, the absolute value of the 
second complex function and the difference between the phase of the wavefront being 
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analyzed and the phase of the second complex function and expressing each of the 
plurality of intensity maps as a sum of the fourth complex function, the fifth complex 
function multiplied by the sine of the known phase delay corresponding to each one of 
the plurality of intensity maps and the sixth complex function multiplied by the cosine 
of the known phase delay corresponding to each one of the plurality of intensity maps. 

Preferably, the step of solving the third function to obtain the amplitude of the 
wavefront being analyzed, the absolute value of the second complex function and the 
difference between the phase of the wavefront being analyzed and the phase of the 
second complex function includes: obtaining two solutions for each of the amplitude of 
the wavefront being analyzed, the absolute value of the second complex function and 
the difference between the phase of the wavefront being analyzed and the phase of the 
second complex function, the two solutions being a higher value solution and a lower 
value solution, combining the two solutions into an enhanced absolute value solution for 
the absolute value of the second complex function, by choosing at each spatial location 
either the higher value solution or the lower value solution of the two solutions in a way 
that the enhanced absolute value solution satisfies the second complex function. 
Preferably, combining the two solutions of the amplitude of the wavefront being 
analyzed into enhanced amplitude solution, by choosing at each spatial location the 
higher value solution or the lower value solution of the two solutions of the amplitude in 
the way that at each location where the higher value solution is chosen for the absolute 
value solution, the higher value solution is chosen for the amplitude solution and at each 
location where the lower value solution is chosen for the absolute value solution, the 
lower value solution is chosen for the amplitude solution, combining the two solutions 
of the difference between the phase of the wavefront being analyzed and the phase of 
the second complex function into an enhanced difference solution, by choosing at each 
spatial location the higher value solution or the lower value solution of the two solutions 
of the difference in the way that at each location where the higher value solution is 
chosen for the absolute value solution, the higher value solution is chosen for the 
difference solution and at each location where the lower value solution is chosen for the 
absolute value solution, the lower value solution is chosen for the difference solution. 

Further in accordance with a preferred embodiment of the present invention the 
spatially uniform, time-varying spatial phase change is applied to a spatially central part 
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of the transformed wavefront. 

Preferably, the transform applied to the wavefront being analyzed is a Fourier 
transform and wherein the step of obtaining a plurality of intensity maps of the plurality 
of phase changed transformed wavefronts includes applying a Fourier transform to the 
pkirality of differently phase changed transformed wavefronts. 

Still further in accordance with a preferred embodiment of the present 
invention the method also includes adding a phase component including relatively high 
frequency components to the wavefront being analyzed prior to applying the transform 
thereto in order to increase the high-frequency content of the transformed wavefront 
prior to the applying the spatially uniform, time-varying spatial phase change to part of 
the transformed wavefront. 

Preferably, the spatially uniform^ time-varying spatial phase change is applied 
to a spatially centered generally circular region of the transformed wavefront and the 
spatially uniform, time-varying spatial phase change is applied to approximately one 
half of the transformed wavefront. 

Additionally in accordance with a preferred embodiment of the present 
invention the transformed wavefront includes a DC region and a non-DC region and the 
spatially uniform, time-varying spatial phase change is applied to at least part of both 
the DC region and the non-DC region. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of differently phase changed transformed wavefronts include a plurality of 
wavefronts whose phase has been changed by employing an at least time varying phase 
change function. Alternatively, the plurality of differently phase changed transformed 
wavefronts include a plurality of wavefronts whose phase has been changed by applying 
an at least time varying phase change function to the wavefront being analyzed. 

Preferably, the at least time varying phase change function is applied to the 
wavefront being analyzed prior to transforming thereof. Alternatively, the at least time 
varying phase change function is applied to the wavefront being analyzed subsequent to 
transforming thereof. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of differently phase changed transformed wavefronts include a plurality of 
wavefronts whose phase has been changed by employing an at least time varying phase 
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change function. 

Additionally or alternatively, the plurality of differently phase changed 
transformed wavefronts include a plurality of wavefronts whose phase has been 
changed by applying an at least time varying phase change function to the wavefront to 
be analyzed. 

Preferably, the at least time varying phase change function is a spatially 
uniform spatial function. 

Additionally in accordance with a preferred embodiment of the present 
invention the transformed wavefront includes a plurality of different wavelength 
components and the plurality of different spatial phase changes are effected by applying 
a phase change to the plurality of different wavelength components of the transformed 
wavefront. 

Preferably, the phase change applied to the plurality of different wavelength 
components of the transformed wavefront is a time-varying spatial phase change. 

Further in accordance with a preferred embodiment of the present invention the 
phase change applied to the plurality of different wavelength components of the 
transformed wavefront is effected by passing the transformed wavefront through an 
object, at least one of whose thickness and refractive index varies spatially. 

Still further in accordance with a preferred embodiment of the present 
invention the phase change applied to the plurality of different wavelength components 
of the transformed wavefront is effected by reflecting the transformed wavefront from a 
spatially varying surface. 

Further in accordance with a preferred embodiment of the present invention the 
phase change applied to the plurality of different wavelength components of the 
transformed wavefront is selected to be different to a predetermined extent for at least 
some of the plurality of different wavelength components. 

Additionally in accordance with a preferred embodiment of the present 
invention the phase change applied to the plurality of different wavelength components 
of the transformed wavefront is identical for at least some of the plurality of different 
wavelength components. 

Further in accordance with a preferred embodiment of the present invention the 
wavefront being analyzed includes a plurality of different wavelength components. 
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Preferably, the plurality of differently phase changed transformed wavefronts are 
obtained by applying a phase change to the plurality of different wavelength 
components of the wavefront being analyzed. 

Preferably, the phase change is applied to the plurality of different wavelength 
components of the wavefront being analyzed prior to transforming thereof. 

Further in accordance with a preferred embodiment of the present invention the 
phase change applied to the plurality of different wavelength components is effected by 
passing the wavefront being analyzed through an object, at least one of whose thickness 
and refractive index varies spatially. 

Further in accordance with a preferred embodiment of the present invention the 
step of obtaining a plurality of intensity maps is performed simultaneously for all of the 
plurality of different wavelength components and obtaining a plurality of intensity maps 
includes dividing the plurality of phase changed transformed wavefronts into separate 
wavelength components. 

Still further in accordance with a preferred embodiment of the present 
invention the step of dividing the plurality of phase changed transformed wavefronts is 
effected by passing the plurality of phase changed transformed wavefronts through a 
dispersion element. 

Additionally in accordance with a preferred embodiment of the present 
invention the phase change applied to the plurality of different wavelength components 
is effected by passing the wavefront being analyzed through an object, at least one of 
whose thickness and refractive index varies spatially, following transforming of the 
wavefront being analyzed. 

Preferably, the phase change which is applied to the plurality of different 
wavelength components is effected by reflecting the wavefront being analyzed from a 
spatially varying surface, following transforming of the wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention the 
phase change applied to the plurality of different wavelength components is selected to 
be different to a predetermined extent for at least some of the plurality of different 
wavelength components. Preferably, the phase change which is applied to the plurality 
of different wavelength components is identical for at least some of the plurality of 
different wavelength components. 



9 



Further in accordance with a preferred embodiment of the present invention the 
phase change applied to the plurality of different wavelength components is effected by 
passing the wavefront being analyzed through a plurality of objects, each characterized 
in that at least one of its thickness and refractive index varies spatially. 

Preferably, the phase change applied to the plurality of different wavelength 
components is effected by passing the wavefront being analyzed through a plurality of 
objects, each characterized in that at least one of its thickness and refractive index varies 
spatially, following transforming of the wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention the 
wavefront being analyzed includes a plurality of different polarization components and 
the plurality of differently phase changed transformed wavefronts are obtained by 
applying a phase change to the plurality of different polarization components of the 
wavefront being analyzed prior to transforming thereof. 

Still further in accordance with a preferred embodiment of the present 
invention the transformed wavefront includes a plurality of different polarization 
components and the plurality of different spatial phase changes are effected by applying 
a phase change to the plurality of different polarization components of the transformed 
wavefront. 

Additionally in accordance with a preferred embodiment of the present 
invention the phase change applied to the plurality of different polarization components 
of the transformed wavefront is different for at least some of the plurality of different 
polarization components. 

Further in accordance with a preferred embodiment of the present invention the 
phase change appHed to the plurality of different polarization components of the 
transformed wavefront is identical for at least some of the plurality of different 
polarization components. 

Additionally in accordance with a preferred embodiment of the present 
invention the step of obtaining a plurality of intensity maps of the plurality of phase 
changed transformed wavefronts includes applying a transform to the plurality of 
differently phase changed transformed wavefronts. 

Preferably, the plurality of phase changed transformed wavefronts are reflected 
from a reflecting surface so that the transform applied to the plurality of differently 
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phase changed transformed wavefronts is identical to the transform applied to the 
wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the wavefront being analyzed is a Fourier transform. 

Still further in accordance with a preferred embodiment of the present 
invention the plurality of intensity maps are obtained by reflecting the plurality of 
differently phase changed transformed wavefronts from a reflecting surface so as to 
transform the plurality of differently phase changed transformed wavefronts. 

Additionally in accordance with a preferred embodiment of the present 
invention the method of obtaining a plurality of intensity maps of the plurality of phase 
changed transformed wavefronts includes applying a transform to the plurality of 
differently phase changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
method of employing the plurality of intensity maps to obtain an output indicating the 
amplitude and phase of the wavefront being analyzed includes expressing the plurality 
of intensity maps as at least one mathematical function of phase and amplitude of the 
wavefront being analyzed and employing the at least one mathematical function to 
obtain an output indicating the phase and amplitude. 

Preferably, the method of employing the plurality of intensity maps to obtain 
an output indicating the amplitude and phase of the wavefront being analyzed includes 
expressing the plurality of intensity maps as at least one mathematical function of phase 
and amplitude of the wavefront being analyzed and of the plurality of different phase 
changes, wherein the phase and amplitude are unknowns and the plurality of different 
phase changes are known and employing the at least one mathematical function to 
obtain an output indicating the phase and amplitude. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of intensity maps includes at least four intensity maps and employing the 
plurality of intensity maps to obtain an output indicating the amplitude and phase of the 
wavefront being analyzed, includes employing a plurality of combinations, each of at 
least three of the plurality of intensity maps, to provide a plurality of indications of the 
amplitude and phase of the wavefront being analyzed. 

Preferably, the method also includes employing the plurality of indications of 
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the amplitude and phase of the wavefront being analyzed to provide an enhanced 
indication of the amplitude and phase of the wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention at 
least some of the plurality of indications of the amplitude and phase are at least second 
order indications of the amplitude and phase of the wavefront being analyzed. 

Further in accordance with a preferred embodiment of the present invention the 
step of obtaining a plurality of differently phase changed transformed wavefronts 
includes applying a transform to the wavefront being analyzed, thereby obtaining a 
transformed wavefront and applying a plurality of different phase and amplitude 
changes to the transformed wavefront, thereby obtaining a plurality of differently phase 
and amplitude changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
plurality of different phase and amplitude changes includes at least three different phase 
and intensity changes, the plurality of different phase and amplitude changes are 
effected by applying at least one of a spatially uniform, time-varying spatial phase 
change and a spatially uniform, time- varying spatial amplitude change to at least part of 
the transformed wavefront, the plurality of intensity maps includes at least three 
intensity maps. Preferably, the step of employing the plurality of intensity maps to 
obtain an output indicating the amplitude and phase of the wavefront being analyzed 
includes: expressing the wavefront being analyzed as a first complex function which has 
an amplitude and phase identical to the amplitude and phase of the wavefront being 
analyzed, expressing the plurality of intensity maps as a function of the first complex 
function and of a spatial function governing at least one of a spatially uniform, 
time-varying spatial phase change and a spatially uniform, time-varying spatial 
amplitude change, defining a second complex function having an absolute value and a 
phase as a convolution of the first complex function and of a Fourier transform of the 
spatial function governing the spatially uniform, time-varying spatial phase change, 
expressing each of the plurality of intensity maps as a third function of: the amplitude of 
the wavefront being analyzed, the absolute value of the second complex function and a 
difference between the phase of the wavefront being analyzed and the phase of the 
second complex function. Preferably, the spatial function governing at least one of a 
spatially uniform, time-varying spatial phase change and a spatially uniform. 
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time-varying spatial amplitude change includes: defining fourth, fifth, sixth and seventh 
complex functions, none of which is a function of any of the plurality of intensity maps 
or of the time-varying spatial phase change. Preferably, each of the fourth, fifth, sixth 
and seventh complex functions being a function of at least one of: the amplitude of the 
wavefront being analyzed, the absolute value of the second complex function and the 
difference between the phase of the wavefront being analyzed and the phase of the 
second complex function, defining an eighth function of a phase delay and of an 
amplitude change, both produced by one of the at least three different phase and 
amplitude changes, corresponding to the at least three intensity maps and expressing 
each of the plurality of intensity maps as a sum of the fourth complex function, the fifth 
complex function multiplied by the absolute value squared of the eighth function, the 
sixth complex function multiplied by the eighth function and the seventh complex 
function multiplied by the complex conjugate of the eighth function, solving the third 
function to obtain the amplitude of the wavefront being analyzed, the absolute value of 
the second complex function and the difference between the phase of the wavefront 
being analyzed and the phase of the second complex function, solving the second 
complex function to obtain the phase of the second complex function and obtaining the 
phase of the wavefront being analyzed by adding the phase of the second complex 
function to the difference between the phase of the wavefront being analyzed and phase 
of the second complex function. 

Further in accordance with a preferred embodiment of the present invention the 
wavefront being analyzed includes at least two wavelength components. Preferably, the 
step of obtaining a plurality of intensity maps also includes dividing the phase changed 
transformed wavefronts according to the at least two wavelength components in order to 
obtain at least two wavelength components of the phase changed transformed 
wavefronts and in order to obtain at least two sets of intensity maps, each set 
corresponding to a different one of the at least two wavelength components of the phase 
changed transformed wavefronts and employing the plurality of intensity maps to obtain 
an output indicating the amplitude and phase of the wavefront being analyzed, obtaining 
an output indicative of the phase of the wavefront being analyzed from each of the at 
least two sets of intensity maps and combining the outputs to provide an enhanced 
indication of phase of the wavefront being analyzed, in which enhanced indication, 
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there is no 2n ambiguity. 

Additionally in accordance with a preferred embodiment of the present 
invention the wavefront being analyzed is an acoustic radiation wavefront. 

Still further in accordance with a preferred embodiment of the present 
invention the wavefront being analyzed includes at least one one-dimensional 
component, the transform applied to the wavefront being analyzed is a one-dimensional 
Fourier transform, performed in a dimension perpendicular to a direction of propagation 
of the wavefront being analyzed, thereby to obtain at least one one-dimensional 
component of the transformed wavefront in the dimension perpendicular to the direction 
of propagation. The plurality of differently phase changed transformed wavefronts are 
obtained by applying the plurality of different phase changes to each of the at least one 
one-dimensional component, thereby obtaining at least one one-dimensional component 
of the plurality of phase changed transformed wavefronts and the plurality of intensity 
maps are employed to obtain an output indicating amplitude and phase of the at least 
one one-dimensional component of the wavefront being analyzed. 

Preferably, the plurality of different phase changes is applied to each of the 
one-dimensional component by providing a relative movement between the wavefront 
being analyzed and an element. Preferably, the element generates spatially varying, 
time-constant phase changes, the relative movement being in an additional dimension 
which is perpendicular both to the direction of propagation and to the dimension 
perpendicular to the direction of propagation. 

Further in accordance with a preferred embodiment of the present invention the 
wavefront being analyzed includes a plurality of different wavelength components, the 
plurality of different phase changes are applied to the plurality of different wavelength 
components of each of the plurality of one-dimensional components of the wavefront 
being analyzed and the step of obtaining a plurality of intensity maps includes dividing 
the plurality of one-dimensional components of the plurality of phase changed 
transformed wavefronts into separate wavelength components. 

Still further in accordance with a preferred embodiment of the present 
invention dividing the plurality of one-dimensional components of the plurality of phase 
changed transformed wavefronts into separate wavelength components is achieved by 
passing the plurality of phase changed transformed wavefronts through a dispersion 
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element. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the wavefront being analyzed includes an additional Fourier 
transform to minimize cross-talk between different one-dimensional components of the 
wavefront being analyzed. 

There is provided in accordance with another preferred embodiment of the 
present invention a method of surface mapping. The method includes obtaining a 
surface mapping wavefront having an amplitude and a phase, by reflecting radiation 
from a surface and analyzing the surface mapping wavefront by: obtaining a plurality of 
differently phase changed transformed wavefronts corresponding to the surface 
mapping wavefront, obtaining a plurality of intensity maps of the plurality of phase 
changed transformed wavefronts and employing the plurality of intensity maps to obtain 
an output indicating the amplitude and phase of the surface mapping wavefront. 

There is further provided in accordance with a preferred embodiment of the 
present invention an apparatus for surface mapping. The apparatus includes a wavefront 
obtainer operating to obtain a surface mapping wavefront having an amplitude and a 
phase, by reflecting radiation from a surface, a wavefront analyzer, analyzing the 
surface mapping wavefront and including a wavefront transformer operating to provide 
a plurality of differently phase changed transformed wavefronts corresponding to the 
surface mapping wavefront, an intensity map generator operating to provide a plurality 
of intensity maps of the plurality of phase changed transformed wavefronts and an 
intensity map utilizer, the plurality of intensity maps provide an output indicating the 
amplitude and phase of the surface mapping wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
radiation reflected from the surface has a narrow band about a given wavelength, 
causing the phase of the surface mapping wavefront to be proportional to geometrical 
variations in the surface, the proportion being an inverse linear function of the 
wavelength. 

Still frirther in accordance with a preferred embodiment of the present 
invention the radiation reflected from the surface has at least two narrow bands, each 
centered about a different wavelength, providing at least two wavelength components in 
the surface mapping wavefront and at least two indications of the phase of the surface 
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mapping wavefront, thereby enabling an enhanced mapping of the surface to be 
obtained by avoiding an ambiguity in the mapping which exceeds the larger of the 
different wavelengths about which the two narrow bands are centered. 

Additionally in accordance with a preferred embodiment of the present 
invention the step of obtaining a plurality of differently phase changed transformed 
wavefronts includes applying a transform to the surface mapping wavefront, thereby to 
obtain a transformed wavefront and applying a plurality of different phase changes, 
including spatial phase changes, to the transformed wavefront, thereby to obtain a 
plurality of differently phase changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the surface mapping wavefront is a Fourier transform, the plurality 
of different phase changes includes at least three different phase changes, effected by 
applying a spatially uniform, time-varying spatial phase change to part of the 
transformed wavefront, the plurality of intensity maps includes at least three intensity 
maps. Preferably, the step of employing the plurality of intensity maps to obtain an 
output indicates the amplitude and phase of the surface mapping wavefront and includes 
expressing the surface mapping wavefront as a first complex function which has an 
amplitude and phase identical to the amplitude and phase of the surface mapping 
wavefront, expressing the plurality of intensity maps as a function of the first complex 
function and of a spatial function governing the spatially imiform, time-varying spatial 
phase change, defining a second complex function having an absolute value and a phase 
as a convolution of the first complex function and of a Fourier transform of the spatial 
function governing the spatially uniform, time-varying spatial phase change. Preferably, 
expressing each of the plurality of intensity maps as a third function of: the amplitude of 
the surface mapping wavefront, the absolute value of the second complex function, a 
difference between the phase of the surface mapping wavefront and the phase of the 
second complex function and a known phase delay produced by one of the at least three 
different phase changes, corresponding to one of the at least three intensity maps, 
solving the third function to obtain the amplitude of the surface mapping wavefront, the 
absolute value of the second complex function and the difference between the phase of 
the surface mapping wavefront and the phase of the second complex function, solving 
the second complex function to obtain the phase of the second complex function and 
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obtaining the phase of the surface mapping wavefront by adding the phase of the second 
complex function to the difference between the phase of the surface mapping wavefront 
and phase of the second complex function. 

Preferably, the surface mapping wavefront includes a plurality of different 
wavelength components. The plurality of differently phase changed transformed 
wavefronts are preferably obtained by: transforming the surface mapping wavefront 
thereby obtaining a transformed wavefront including a plurality of different wavelength 
components and applying a phase change to the plurality of different wavelength 
components of the transformed wavefront by passing the transformed wavefront 
through an object, at least one of whose thickness and refractive index varies spatially. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention a method of inspecting an object. The method includes obtaining 
an object inspection wavefront which has an amplitude and a phase, by transmitting 
radiation through the object and analyzing the object inspection wavefront by: obtaining 
a plurality of differently phase changed transformed wavefronts corresponding to the 
object inspection wavefront, obtaining a plurality of intensity maps of the plurality of 
phase changed transformed wavefronts and employing the plurality of intensity maps to 
obtain an output indicating the amplitude and phase of the object inspection wavefront. 

There is farther provided in accordance with a preferred embodiment of the 
present invention an apparatus for inspecting an object. The apparatus includes a 
wavefront obtainer operating to obtain an object inspection wavefront which has an 
amplitude and a phase, by transmitting radiation through the object, a wavefront 
analyzer, analyzing the object inspection wavefront and including a wavefront 
transformer operating to provide a plurality of differently phase changed transformed 
wavefronts corresponding to the object inspection wavefront, an intensity map generator 
operating to provide a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts and an intensity map utilizer, employing the plurality of 
intensity maps to provide an output indicating the amplitude and phase of the object 
inspection wavefront. 

Preferably, when the object is substantially uniform in material and other 
optical properties, the phase of the object inspection wavefront is proportional to the 
object thickness. 
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Additionally, when the object is substantially uniform in thickness, the phase 
of the object inspection wavefront is proportional to optical properties of the object. 

Further in accordance with a preferred embodiment of the present invention the 
radiation has at least two narrow bands, each centered about a different wavelength, 
providing at least two wavelength components in the object inspection wavefront and at 
least two indications of the phase of the object inspection wavefront, thereby enabling 
an enhanced mapping of thickness of the object to be inspected by avoiding an 
ambiguity in the mapping which exceeds the larger of the different wavelengths about 
which the two narrow bands are centered. 

Still further in accordance with a preferred embodiment of the present 
invention the method of obtaining a plurality of differently phase changed transformed 
wavefronts includes applying a transform to the object inspection wavefront, thereby 
obtaining a transformed wavefront and applying a plurality of different phase changes, 
including spatial phase changes, to the transformed wavefront, thereby obtaining a 
plurality of differently phase changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the object inspection wavefront is a Fourier transform, the plurality 
of different phase changes includes at least three different phase changes, effected by 
applying a spatially uniform, time-varying spatial phase change to part of the 
transformed wavefront. Preferably, the plurality of intensity maps includes at least three 
intensity maps and employing the plurality of intensity maps to obtain an output 
indicating the amplitude and phase of the object inspection wavefront and includes: 
expressing the object inspection wavefront as a first complex fimction which has an 
amplitude and phase identical to the amplitude and phase of the object inspection 
wavefront, expressing the plurality of intensity maps as a frmction of the first complex 
function and of a spatial frinction governing the spatially uniform, time-varying spatial 
phase change, defining a second complex function having an absolute value and a phase 
as a convolution of the first complex function and of a Fourier transform of the spatial 
function governing the spatially uniform, time-varying spatial phase change, expressing 
each of the plurality of intensity maps as a third frmction of: the amplitude of the object 
inspection wavefront, the absolute value of the second complex fimction, a difference 
between the phase of the object inspection wavefront and the phase of the second 
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complex function and a known phase delay produced by one of the at least three 
different phase changes, corresponding to one of the at least three intensity maps, 
solving the third function to obtain the amplitude of the object inspection waveftont, the 
absolute value of the second complex function and the difference between the phase of 
the object inspection wave front and the phase of the second complex function, solving 
the second complex function to obtain the phase of the second complex function and 
obtaining the phase of the object inspection wavefront by adding the phase of the 
second complex function to the difference between the phase of the object inspection 
wavefront and phase of the second complex function. 

Still further in accordance with a preferred embodiment of the present 
invention the object inspection wavefront includes a plurality of different wavelength 
components. The plurality of differently phase changed transformed wavefronts are 
preferably obtained by: transforming the object inspection wavefront thereby obtaining 
a transformed wavefront including a plurality of different wavelength components and 
applying a phase change to the plurality of different wavelength components of the 
transformed wavefront by reflecting the transformed wavefront from a spatially varying 
surface. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention a method of spectral analysis. The method includes obtaining a 
spectral analysis wavefront having an amplitude and a phase, by causing radiation to 
impinge on an object, analyzing the spectral analysis wavefront by: obtaining a plurality 
of differently phase changed transformed wavefronts corresponding to the spectral 
analysis wavefront which has an amplitude and a phase, obtaining a plurality of 
intensity maps of the plurality of phase changed transformed wavefronts and employing 
the plurality of intensity maps to obtain an output indicating the amplitude and phase of 
the spectral analysis wavefront and employing the output indicating the amplitude and 
phase to obtain an output indicating spectral content of the radiation. 

There is provided in accordance with a further preferred embodiment of the 
present invention an apparatus for spectral analysis. The apparatus includes a wavefront 
obtainer operating to obtain a spectral analysis wavefront having an amplitude and a 
phase, by causing radiation to impinge on an object, a wavefront analyzer, analyzing the 
spectral analysis wavefront, including a wavefront transformer operating to provide a 
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plurality of differently phase changed transformed wavefronts corresponding to the 
spectral analysis wavefront which has an amplitude and a phase, an intensity map 
generator operating to provide a plurality of intensity maps of the plurality of phase 
changed transformed wavefronts, an intensity map utilizer, employing the plurality of 
intensity maps to provide an output indicating the amplitude and phase of the spectral 
analysis wavefront and a phase and amplitude utilizer, employing the output indicating 
the amplitude and phase to obtain an output indicating spectral content of the radiation. 

Further in accordance with a preferred embodiment of the present invention 
and wherein obtaining the spectral analysis wavefront is effected by reflecting the 
radiation from the object. 

Still further in accordance with a preferred embodiment of the present 
invention and wherein obtaining the spectral analysis wavefront is effected by 
transmitting the radiation through the object. 

Additionally in accordance with a preferred embodiment of the present 
invention the radiation is substantially of a single wavelength, the phase of the spectral 
analysis wavefront is inversely proportional to the single wavelength, and is related to at 
least one of a surface characteristic and thickness of the impinged object. 

Still further in accordance with a preferred embodiment of the present 
invention the step of employing the plurality of intensity maps to obtain an output 
indicating the amplitude and phase of the spectral analysis wavefront includes: 
expressing the plurality of intensity maps as at least one mathematical function of phase 
and amplitude of the spectral analysis wavefront and of the plurality of different phase 
changes, wherein at least the phase is unknown and a function generating the plurality 
of phase changed transformed wavefronts is known and employing the at least one 
mathematical function to obtain an output indicating at least the phase. 

Additionally in accordance with a preferred embodiment of the present 
invention the step of obtaining a plurality of differently phase changed transformed 
wavefronts includes applying a transform to the spectral analysis wavefront, thereby 
obtaining a transformed wavefront and applying a plurality of different phase changes, 
including spatial phase changes, to the transformed wavefront, thereby obtaining a 
plurality of differently phase changed transformed wavefronts. 
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Further in accordance with a preferred embodiment of the present invention the 
transform applied to the spectral analysis wavefront is a Fourier transform, the plurality 
of different phase changes includes at least three different phase changes, effected by 
applying a spatially uniform, time-varying spatial phase change to part of the 
transformed v^avefront. Preferably, the plurality of intensity maps includes at least three 
intensity maps and employing the plurality of intensity maps to obtain an output 
indicating the amplitude and phase of the spectral analysis wavefront includes: 
expressing the spectral analysis wavefront as a first complex function which has an 
amplitude and phase identical to the amplitude and phase of the spectral analysis 
wavefront, expressing the plurality of intensity maps as a function of the first complex 
function and of a spatial function governing the spatially uniform, time-varying spatial 
phase change, defining a second complex function having an absolute value and a phase 
as a convolution of the first complex function and of a Fourier transform of the spatial 
function governing the spatially uniform, time-varying spatial phase change, expressing 
each of the plurality of intensity maps as a third function of: the amplitude of the 
spectral analysis wavefront, the absolute value of the second complex function, a 
difference between the phase of the spectral analysis wavefront and the phase of the 
second complex function and a known phase delay produced by one of the at least three 
different phase changes, corresponding to one of the at least three intensity maps, 
solving the third function to obtain the amplitude of the spectral analysis wavefront, the 
absolute value of the second complex function and the difference between the phase of 
the spectral analysis wavefront and the phase of the second complex function, solving 
the second complex function to obtain the phase of the second complex function and 
obtaining the phase of the spectral analysis wavefront by adding the phase of the second 
complex function to the difference between the phase of the spectral analysis wavefront 
and phase of the second complex function. 

Further in accordance with a preferred embodiment of the present invention the 
spectral analysis wavefront includes a plurality of different wavelength components and 
the plurality of differently phase changed transformed wavefronts are obtained by 
applying a phase change to the plurality of different wavelength components of the 
spectral analysis wavefront. 

There is further provided in accordance with a preferred embodiment of the 
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present invention a method of phase change analysis. The method includes obtaining a 
phase change analysis wavefront which has an amplitude and a phase, applying a 
transform to the phase change analysis wavefront thereby to obtain a transformed 
wavefront, applying a plurality of different phase changes to the transformed wavefront, 
thereby to obtain a plurality of differently phase changed transformed wavefronts, 
obtaining a plurality of intensity maps of the plurality of phase changed transformed 
wavefronts and employing the plurality of intensity maps to obtain an output indication 
of differences between the plurality of different phase changes applied to the 
transformed phase change analysis wavefront. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention an apparatus for phase change analysis. The apparatus includes a 
wavefront obtainer, operating to obtain a phase change analysis wavefront which has an 
amplitude and a phase, a transform applier, applying a transform to the phase change 
analysis wavefront thereby to obtain a transformed wavefront, a phase change applier, 
applying at least one phase change to the transformed wavefront, thereby to obtain at 
least one phase changed transformed wavefront, an intensity map generator operating to 
provide at least one intensity map of the phase changed transformed wavefront and an 
intensity map utilizer, employing the plurality of intensity maps to provide an output 
indication of differences between the plurality of different phase changes applied to the 
transformed phase change analysis wavefront. 

Typically, when lateral shifts appear in the plurality of different phase changes, 
corresponding changes appear in the plurality of intensity maps and the step of 
employing the plurality of intensity maps results in obtaining an indication of the lateral 
shifts. 

Still further in accordance with a preferred embodiment of the present 
invention the step of employing the plurality of intensity maps to obtain an output 
indication of differences between the plurality of different phase changes applied to the 
transformed phase change analysis wavefront includes: expressing the plurality of 
intensity maps as at least one mathematical function of phase and amplitude of the 
phase change analysis wavefront and of the plurality of different phase changes, where 
at least the phase and amplitude are known and the plurality of different phase changes 
are unknown and employing the mathematical function to obtain an output indicating 
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the differences between the plurality of different phase changes. 

There is further provided in accordance with yet a further preferred 
embodiment of the present invention a method of phase change analysis. The method 
includes obtaining a phase change analysis wavefront which has an amplitude and a 
phase, applying a transform to the phase change analysis wavefront thereby to obtain a 
transformed wavefront, applying at least one phase change to the transformed 
wavefront, thereby to obtain at least one phase changed transformed wavefront, 
obtaining at least one intensity map of the at least one phase changed transformed 
wavefront and employing the intensity map to obtain an output indication of the at least 
one phase change applied to the transformed phase change analysis wavefront. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention an apparatus for phase change analysis. The apparatus includes a 
wavefront obtainer, operating to obtain a phase change analysis wavefront which has an 
amplitude and a phase, a transform applier, applying a transform to the phase change 
analysis wavefront thereby to obtain a transformed wavefront, a phase change applier, 
applying at least one phase change to the transformed wavefront, thereby to obtain at 
least one phase changed transformed wavefront, an intensity map generator operating to 
provide at least one intensity map of the phase changed transformed wavefront and an 
intensity map utilizer, employing the intensity map to provide an output indication of 
the phase change applied to the transformed phase change analysis wavefront. 

Preferably, the phase change is a phase delay, having a value selected from a 
plurality of pre-determined values, and the output indication of the phase change 
includes the value of the phase delay. 

There is also provided in accordance with a preferred embodiment of the 
present invention a method of stored data retrieval. The method includes obtaining a 
stored data retrieval wavefront which has an amplitude and a phase, by reflecting 
radiation from the media in which information is encoded, by selecting the height of the 
media at each of a multiplicity of different locations on the media. Preferably, analyzing 
the stored data retrieval wavefront by: obtaining a plurality of differently phase changed 
transformed wavefronts corresponding to the stored data retrieval wavefront, obtaining 
a plurality of intensity maps of the plurality of phase changed transformed wavefronts 
and employing the plurality of intensity maps to obtain an indication of the amplitude 
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and phase of the stored data retrieval wavefront and employing the indication of the 
amplitude and phase to obtain the information. 

There is further provided in accordance with yet another preferred embodiment 
of the present invention an apparatus for stored data retrieval. The apparatus includes a 
wavefront obtainer operating to obtain a stored data retrieval wavefront which has an 
amplitude and a phase, by reflecting radiation from the media in which information is 
encoded by selecting the height of the media at each of a multiplicity of different 
locations on the media, a wavefront analyzer, analyzing the stored data retrieval 
wavefront and including a wavefront transformer operating to provide a plurality of 
differently phase changed transformed wavefronts corresponding to the stored data 
retrieval wavefront, an intensity map generator operating to obtain a plurality of 
intensity maps of the plurality of phase changed transformed wavefronts and an 
intensity map utilizer, employing the plurality of intensity maps to provide an indication 
of the amplitude and phase of the stored data retrieval wavefront and a phase and 
amplitude utilizer, employing the indication of the amplitude and phase to provide the 
information. 

Preferably, the step of obtaining a plurality of differently phase changed 
transformed wavefronts includes: applying a transform to the stored data retrieval 
wavefront thereby to obtain a transformed wavefront and applying a plurality of 
different phase changes to the transformed wavefront, thereby to obtain a plurality of 
differently phase changed transformed wavefronts. 

Further in accordance with a preferred embodiment of the present invention the 
transform applied to the stored data retrieval wavefront is a Fourier transform, the 
plurality of different phase changes includes at least three different phase changes, 
effected by applying a spatially uniform, time-varying spatial phase change to part of 
the transformed wavefront, the plurality of intensity maps includes at least three 
intensity maps and employing the plurality of intensity maps to obtain an output 
indicating the amplitude and phase of the stored data retrieval wavefront includes: 
expressing the stored data retrieval wavefront as a first complex function which has an 
amplitude and phase identical to the amplitude and phase of the stored data retrieval 
wavefront, expressing the plurality of intensity maps as a function of the first complex 
function and of a spatial function governing the spatially uniform, time-varying spatial 
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phase change, defining a second complex function having an absolute value and a phase 
as a convolution of the first complex function and of a Fourier transform of the spatial 
function governing the spatially uniform, time-varying spatial phase change, expressing 
each of the plurality of intensity maps as a third function of: the amplitude of the stored 
data retrieval wavefront, the absolute value of the second complex function, a difference 
between the phase of the stored data retrieval wavefront and the phase of the second 
complex function and a known phase delay produced by one of the at least three 
different phase changes, corresponding to one of the at least three intensity maps, 
solving the third function to obtain the amplitude of the stored data retrieval wavefront, 
the absolute value of the second complex function and the difference between the phase 
of the stored data retrieval wavefront and the phase of the second complex function, 
solving the second complex function to obtain the phase of the second complex function 
and obtaining the phase of the stored data retrieval wavefront by adding the phase of the 
second complex function to the difference between the phase of the stored data retrieval 
wavefront and phase of the second complex function. 

Still further in accordance with a preferred embodiment of the present 
invention the stored data retrieval wavefront includes at least one one-dimensional 
component, the transform applied to the data retrieval wavefront is a one-dimensional 
Fourier transform, performed in a dimension perpendicular to a direction of propagation 
of the data retrieval wavefront, thereby to obtain at least one one-dimensional 
component of the transformed wavefront in the dimension perpendicular to the direction 
of propagation, the plurality of differently phase changed transformed wavefironts are 
obtained by applying the plurality of different phase changes to each of the 
one-dimensional component, thereby obtaining at least one one-dimensional component 
of the plurality of phase changed transformed wavefronts and the plurality of intensity 
maps are employed to obtain an output indicating amplitude and phase of the 
one-dimensional component of the data retrieval wavefront. 

Preferably, the plurality of different phase changes is applied to each of the at 
least one one-dimensional component by providing a relative movement between the 
media and a component generating spatially varying, time-constant phase changes, the 
relative movement being in a dimension perpendicular to the direction of propagation 
and to the dimension perpendicular to the direction of propagation. 
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Additionally in accordance with a preferred embodiment of the present 
invention the information is encoded on the media whereby: an intensity value is 
realized by reflection of light from each location on the media to lie within a 
predetermined range of values, the range corresponding an element of the information 
stored at the location and by employing the plurality of intensity maps, multiple 
intensity values are realized for each location, providing multiple elements of 
information for each location on the media. 

Preferably, the plurality of differently phase changed transformed wavefronts 
include a plurality of wavefronts whose phase has been changed by applying an at least 
time varying phase change function to the stored data retrieval wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
stored data retrieval wavefront includes a plurality of different wavelength components 
and the plurality of differently phase changed transformed wavefronts are obtained by 
applying at least one phase change to the plurality of different wavelength components 
of the stored data retrieval wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
radiation which is reflected from the media includes a plurahty of different wavelength 
components, resulting in the stored data retrieval wavefront including a plurality of 
different wavelength components and the plurality of differently phase changed 
transformed wavefronts are obtained by applying a phase change to the plurality of 
different wavelength components of the stored data retrieval wavefront. 

Still further in accordance v^th a preferred embodiment of the present 
invention the information encoded by selecting the height of the media at each of a 
multiplicity of different locations on the media is also encoded by selecting the 
reflectivity of the media at each of a plurality of different locations on the media and 
employing the indication of the amplitude and phase to obtain the information includes 
employing the indication of the phase to obtain the information encoded by selecting the 
height of the media and employing the indication of the amplitude to obtain the 
information encoded by selecting the reflectivity of the media. 

There is provided in accordance with another preferred embodiment of the 
present invention a method of 3 -dimensional imaging. The method includes obtaining a 
3 -dimensional imaging wavefront, which has an amplitude and a phase, by reflecting 
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radiation from an object to be viewed and analyzing the 3 -dimensional imaging 
wavefront by: obtaining a plurality of differently phase changed transformed wavefronts 
corresponding to the 3-dimensional imaging wavefront, obtaining a plurality of intensity 
maps of the plurality of differently phase changed transformed wavefronts and 
employing the plurality of intensity maps to obtain an output indicating the amplitude 
and phase of the 3-dimensional imaging wavefront. 

There is further provided in accordance with a preferred embodiment of the 
present invention an apparatus for 3-dimensional imaging. The apparatus includes a 
wavefront obtainer operating to obtain a 3-dimensional imaging wavefront, which has 
an amplitude and a phase, by reflecting radiation from an object to be viewed, a 
wavefront analyzer, analyzing the 3-dimensional imaging wavefront including a 
wavefront transformer operative to provide a plurality of differently phase changed 
transformed wavefronts corresponding to the 3-dimensional imaging wavefront, an 
intensity map generator operative to provide a plurality of intensity maps of the plurality 
of differently phase changed transformed wavefronts and an intensity map utilizer, 
employing the plurality of intensity maps to provide an output indicating the amplitude 
and phase of the 3-dimensional imaging wavefront. 

Further in accordance with a preferred embodiment of the present invention the 
radiation reflected from the object has a narrow band about a given wavelength, causing 
the phase of the 3-dimensional imaging wavefront to be proportional to geometrical 
variations in the object, the proportion being an inverse linear function of the 
wavelength. 

Additionally in accordance with a preferred embodiment of the present 
invention the step of obtaining a plurality of differently phase changed transformed 
wavefronts includes applying a transform to the 3-dimensional imaging wavefront, 
thereby to obtain a transformed wavefront and applying a plurality of different phase 
changes, including spatial phase changes, to the transformed wavefront, thereby to 
obtain a plurality of differently phase changed transformed wavefronts. 

Still further in accordance with a preferred embodiment of the present 
invention the 3-dimensional imaging wavefront includes a plurality of different 
wavelength components and the plurality of differently phase changed transformed 
wavefronts are obtained by: transforming the 3-dimensional imaging wavefront, thereby 
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obtaining a transformed wavefront including a plurality of different wavelength 
components and applying phase changes to the plurality of different wavelength 
components of the transformed wavefront by passing the transformed wavefront 
through an object, at least one of whose thickness and refractive index varies spatially. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention a method of wavefront analysis. The method includes obtaining a 
plurality of differently phase changed transformed wavefronts corresponding to a 
wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of 
phase changed transformed wavefronts and employing the plurality of intensity maps to 
obtain an output indicating at least the phase of the wavefront being analyzed by 
combining the plurality of intensity maps into a second plurality of combined intensity 
maps, the second plurality being less than the first plurality, obtaining at least an output 
indicative of the phase of the wavefront being analyzed from each of the second 
plurality of combined intensity maps and combining the outputs to provide at least an 
enhanced indication of phase of the wavefront being analyzed. 

There is also provided in accordance with yet another preferred embodiment of 
the present invention an apparatus wavefront analysis. The apparatus includes a 
wavefront transformer operating to provide a plurality of differently phase changed 
transformed wavefronts corresponding to a wavefront being analyzed, an intensity map 
generator operating to obtain a plurality of intensity maps of the plurality of phase 
changed transformed wavefronts and an intensity map utilizer, employing the plurality 
of intensity maps to obtain an output indicating at least amplitude of the wavefront 
being analyzed and including an intensity combiner operating to combine the plurality 
of intensity maps into a second plurality of combined intensity maps, the second 
plurality being less than the first plurality, an indication provider operating to provide at 
least an output indicative of the amplitude of the wavefront being analyzed from each of 
the second plurality of combined intensity maps and an enhanced indication provider, 
combining the outputs to provide at least an enhanced indication of amplitude of the 
wavefront being analyzed. 

There is provided in accordance with a further preferred embodiment of the 
present invention a method of wavefront analysis. The method includes obtaining a 
plurality of differently phase changed transformed wavefronts corresponding to a 
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wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of 
phase changed transformed wavefront and employing the plurality of intensity maps to 
obtain an output indicating at least amplitude of the wavefront being analyzed by 
combining the plurality of intensity maps into a second plurality of combined intensity 
maps, the second plurality being less than the first plurality, obtaining at least an output 
indicative of the amplitude of the wavefront being analyzed from each of the second 
plurality of combined intensity maps and combining the outputs to provide at least an 
enhanced indication of amplitude of the wavefront being analyzed. 

There is provided in accordance with a preferred embodiment of the present 
invention an apparatus for wavefront analysis. The apparatus includes a wavefront 
transformer operating to provide a plurality of differently phase changed transformed 
wavefronts corresponding to a wavefront being analyzed, an intensity map generator 
operating to provide a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts and an intensity map utilizer, employing the plurality of 
intensity maps to provide an output indicating at least the phase of the wavefront being 
analyzed. Preferably, the apparatus also includes an intensity map expresser, expressing 
the plurality of intensity maps as a fimction of: amplitude of the wavefront being 
analyzed, phase of the wavefront being analyzed and a phase change function 
characterizing the plurality of differently phase changed transformed wavefronts, a 
complex function definer, defining a complex function of: the amplitude of the 
wavefront being analyzed, the phase of the wavefront being analyzed and the phase 
change function characterizing the plurality of differently phase changed transformed 
wavefronts, the complex function being characterized in that the intensity at each 
location in the plurality of intensity maps is a function predominantly of a value of the 
complex fimction at the location and of the amplitude and the phase of the wavefront 
being analyzed at the location. The apparatus also typically, includes complex fiinction 
expresser, expressing the complex function as a fimction of the plurality of intensity 
maps and a phase obtainer, obtaining values for the phase by employing the complex 
function expressed as a fiinction of the plurality of intensity maps. 

There is also provided in accordance with another preferred embodiment of the 
present invention a method of wavefront analysis. The method includes obtaining a 
plurality of differently phase changed transformed wavefronts corresponding to a 



29 



wavefront being analyzed, obtaining a plurality of intensity maps of the plurality of 
phase changed transformed wavefronts and employing the plurality of intensity maps to 
provide an output indicating at least the phase of the wavefront being analyzed by: 
expressing the plurality of intensity maps as a function of: amplitude of the wavefront 
being analyzed, phase of the wavefront being analyzed and a phase change function 
characterizing the plurality of differently phase changed transformed wavefronts. 
Additionally, defining a complex function of: the amplitude of the wavefront being 
analyzed, the phase of the wavefront being analyzed and the phase change function 
characterizing the plurality of differently phase changed transformed wavefronts, the 
complex function being characterized in that the intensity at each location in the 
plurality of intensity maps is a function predominantly of a value of the complex 
function at the location and of the amplitude and the phase of the wavefront being 
analyzed at the location, expressing the complex function as a function of the plurality 
of intensity maps and obtaining values for the phase by employing the complex function 
expressed as a function of the plurality of intensity maps. 

There is further provided in accordance with yet a further preferred 
embodiment of the present invention a method of wavefront analysis. The method 
includes applying a Fourier transform to a wavefront being analyzed which has an 
amplitude and a phase, thereby obtaining a transformed wavefront, applying a spatially 
uniform time-varying spatial phase change to part of the transformed wavefront, thereby 
to obtain at least three differently phase changed transformed wavefronts, applying a 
second Fourier transform to obtain at least three intensity maps of the at least three 
phase changed transformed wavefronts and employing the at least three intensity maps 
to obtain an output indicating at least one of the phase and the amplitude of the 
wavefront being analyzed by: expressing the wavefront being analyzed as a first 
complex function which has an amplitude and phase identical to the amplitude and 
phase of the wavefront being analyzed, expressing the plurality of intensity maps as a 
function of the first complex function and of a spatial function governing the spatially 
uniform, time-varying spatial phase change, defining a second complex function having 
an absolute value and a phase as a convolution of the first complex function and of a 
Fourier transform of the spatial function governing the spatially uniform, time-varying 
spatial phase change, expressing each of the plurality of intensity maps as a third 
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function of: the amplitude of the wavefront being analyzed, the absolute value of the 
second complex function, a difference between the phase of the wavefront being 
analyzed and the phase of the second complex function and a known phase delay 
produced by one of the at least three different phase changes, which each correspond to 
one of the at least three intensity maps, solving the third function to obtain the 
amplitude of the wavefront being analyzed, the absolute value of the second complex 
function and the difference between the phase of the wavefront being analyzed and the 
phase of the second complex function, solving the second complex function to obtain 
the phase of the second complex function and obtaining the phase of the wavefront 
being analyzed by adding the phase of the second complex function to the difference 
between the phase of the wavefront being analyzed and phase of the second complex 
function. 

There is further provided in accordance with yet a further preferred 
embodiment of the present invention an apparatus for wavefront analysis. The apparatus 
includes a first transform applier, applying a Fourier transform to a wavefront being 
analyzed which has an amplitude and a phase thereby to obtain a transformed 
wavefront, a phase change applier, applying a spatially uniform time-varying spatial 
phase change to part of the transformed wavefront, thereby obtaining at least three 
differently phase changed transformed wavefronts, a second transform applier, applying 
a second Fourier transform to the at least three differently phase changed transformed 
wavefronts, thereby obtaining at least three intensity maps. The apparatus also typically 
includes an intensity map utilizer, employing the at least three intensity maps to provide 
an output indicating the phase and the amplitude of the wavefront being analyzed and a 
wavefront expresser, expressing the wavefront being analyzed as a first complex 
function which has an amplitude and phase identical to the amplitude and phase of the 
wavefront being analyzed, a first intensity map expresser, expressing the plurality of 
intensity maps as a function of the first complex function and of a spatial function 
governing the spatially uniform, time-varying spatial phase change. Preferably, the 
apparatus also includes a complex function definer, defining a second complex function 
having an absolute value and a phase as a convolution of the first complex function and 
of a Fourier transform of the spatial function governing the spatially uniform, 
time- varying spatial phase change, a second intensity map expresser, expressing each of 
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the plurality of intensity maps as a third function of; the amplitude of the wavefront 
being analyzed, the absolute value of the second complex function, a difference between 
the phase of the wavefront being analyzed and the phase of the second complex function 
and a knoAvn phase delay produced by one of the at least three different phase changes, 
which each correspond to one of the at least three intensity maps. The apparatus further 
typically includes a first function solver, solving the third function to obtain the 
amplitude of the wavefront being analyzed, the absolute value of the second complex 
function and the difference between the phase of the wavefront being analyzed and the 
phase of the second complex function, a second function solver, solving the second 
complex function to obtain the phase of the second complex function and a phase 
obtainer, obtaining the phase of the wavefront being analyzed by adding the phase of 
the second complex function to the difference between the phase of the wavefront being 
analyzed and the phase of the second complex function. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention will be understood and appreciated more fully from the 
following detailed description, taken in conjunction with the drawings in which: 

Fig. lA is a simplified partially schematic, partially pictorial illustration of 
wavefront analysis functionality operative in accordance with a preferred embodiment 
of the present invention; 

Fig, IB is a simplified partially schematic, partially block diagram illustration 
of a wavefront analysis system suitable for carrying out the functionality of Fig. lA in 
accordance with a preferred embodiment of the present invention; 

Fig, 2 is a simplified functional block diagram illustration of the functionality 
of Fig. 1 A where time-varying phase changes are applied to a transformed wavefront; 

Fig. 3 is a simplified functional block diagram illustration of the functionality 
of Fig. lA where time-varying phase changes are applied to a wavefront prior to 
transforming thereof; 

Fig. 4 is a simplified functional block diagram illustration of the functionality 
of Fig. 2 where time-varying, non-spatially varying spatial phase changes are applied to 
a transformed wavefront; 

Fig. 5 is a simplified functional block diagram illustration of the functionality 
of Fig. 3 where time-varying, non-spatially varying spatial phase changes are applied to 
a wavefront prior to transforming thereof; 

Fig. 6 is a simplified functional block diagram illustration of the functionality 
of Fig. lA where phase changes are applied to a plurality of different wavelength 
components of a transformed wavefront; 

Fig. 7 is a simplified functional block diagram illustration of the functionality 
of Fig. lA where phase changes are applied to a plurality of different wavelength 
components of a wavefront prior to transforming thereof; 

Fig. 8 is a simplified functional block diagram illustration of the functionality 
of Fig. lA where phase changes are applied to a plurality of different polarization 
components of a transformed wavefront; 

Fig. 9 is a simplified functional block diagram illustration of the functionality 
of Fig. lA where phase changes are applied to a plurality of different polarization 
components of a wavefront prior to transforming thereof; 
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Fig. lOA is a simplified functional block diagram illustration of the 
functionality of Fig. lA where a wavefront being analyzed comprises at least one 
one-dimensional component; 

Fig. 1 OB is a simplified partially schematic, partially pictorial illustration of a 
wavefront analysis system suitable for carrying out the functionality of Fig. lOA in 
accordance with a preferred embodiment of the present invention; 

Fig. 1 1 is a simplified functional block diagram illustration of the functionality 
of Fig. I A where an additional transform is applied following the application of spatial 
phase changes; 

Fig. 12 is a simplified functional block diagram illustration of the functionality 
of Fig. lA, wherein intensity maps are employed to provide information about a 
wavefront being analyzed, such as indications of amplitude and phase of the wavefront; 

Fig. 13 is a simplified functional block diagram illustration of part of the 
functionality of Fig. lA, wherein the transform applied to the wavefront being analyzed 
is a Fourier transform, wherein at least three different spatial phase changes are applied 
to a transformed wavefront. and wherein at least three intensity maps are employed to 
obtain indications of at least the phase of a wavefront; 

Fig. 14 is a simplified partially schematic, partially pictorial illustration of part 
of one preferred embodiment of a wavefront analysis system of the type shown in Fig. 
IB; 

Fig. 15 is a simplified partially schematic, partially pictorial illustration of a 
system for surface mapping employing the functionality and structure of Figs. 1 A and 
IB; 

Fig. 16 is a simplified partially schematic, partially pictorial illustration of a 
system for object inspection employing the functionality and structure of Figs. lA and 
IB; 

Fig, 17 is a simplified partially schematic, partially pictorial illustration of a 
system for spectral analysis employing the fionctionality and structure of Figs. 1 A and 
IB; 

Fig. 18 is a simplified partially schematic, partially pictorial illustration of a 
system for phase-change analysis employing the functionality and structure of Figs. lA 
and IB; 
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Fig, 19 is a simplified partially schematic, partially pictorial illustration of a 
system for stored data retrieval employing the functionality and structure of Figs. 1 A 
and IB; 

Fig. 20 is a simplified partially schematic, partially pictorial illustration of a 
system for 3-dimensional imaging employing the functionality and structure of Figs. lA 
and IB; 

Fig. 21 A is a simplified partially schematic, partially pictorial illustration of 
wavefront analysis functionality operative in accordance with another preferred 
embodiment of the present invention; 

Fig. 21 B is a simplified partially schematic, partially block diagram illustration 
of a wavefront analysis system suitable for carrying out the functionalit)^ of Fig. 21 A in 
accordance with another preferred embodiment of the present invention; and 

Fig. 22 is a simplified partially schematic, partially pictorial illustration of a 
system for surface mapping employing the functionality and structure of Figs. 21 A and 
21B. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Reference is now made to Fig. lA, which is a simplified partially schematic, 
partially pictorial illustration of wavefront analysis functionality operative in 
accordance with a preferred embodiment of the present invention. The functionality of 
Fig. 1 A can be summarized as including the following sub-functionalities: 

A. obtaining a plurality of differently phase changed transformed wavefronts 
corresponding to a wavefront being analyzed, which has an amplitude and a phase; 

B. obtaining a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts; and 

C. employing the plurality of intensity maps to obtain an output indicating at least 
one and possibly both of the phase and the amplitude of the wavefront being analyzed. 

As seen in Fig. lA, the first sub-functionality, designated "A" may be realized 
by the following functionalities: 

A wavefront, which may be represented by a plurality of point sources of light, 
is generally designated by reference numeral 100. Wavefront 100 has a phase 
characteristic which is typically spatially non-uniform, shown as a solid line and 
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indicated generally by reference numeral 102. Wavefront 100 also has an amplitude 
characteristic which is also typically spatially non-uniform, shown as a dashed line and 
indicated generally by reference numeral 103. Such a wavefront may be obtained in a 
conventional manner by receiving Ught from any object, such as by reading an optical 
disk, for example a DVD or compact disk 104. 

A principal purpose of the present invention is to measure the phase 
characteristic, such as that indicated by reference numeral 102, which is not readily 
measured. Another purpose of the present invention is to measure the amplitude 
characteristic, such as that indicated by reference numeral 103 in an enhanced manner, 
A further purpose of the present invention is to measure both the phase characteristic 
102 and the amplitude characteristic 103. While there exist various techniques for 
carrying out such measurements, the present invention provides a methodology which is 
believed to be superior to those presently known, inter aha due to its relative 
insensitivity to noise. 

A transform, indicated here symbolically by reference numeral 106, is applied 
to the wavefront being analyzed 100, thereby to obtain a transformed wavefront. A 
preferred transform is a Fourier transform. The resulting transformed wavefront is 
symbolically indicated by reference numeral 108. 

A plurality of different phase changes, preferably spatial phase changes, 
represented by optical path delays 110, 112 and 114 are apphed to the transformed 
wavefront 108, thereby to obtain a plurality of differently phase changed transformed 
wavefronts, represented by reference numerals 120, 122 and 124 respectively. It is 
appreciated that the illustrated difference between the individual ones of the plurality of 
differently phase changed transformed wavefronts is that portions of the transformed 
wavefront are delayed differently relative to the remainder thereof. The difference in the 
phase changes, which are applied to the transformed wavefront 108, is represented in 
Fig. 1 A by the change in thickness of the optical path delays 110, 1 12 and 1 14. 

As seen in Fig. lA, the second sub-fimctionality, designated "B", may be 
realized by applying a transform, preferably a Fourier transform, to the plurality of 
differently phase changed transformed wavefronts. Alternatively, the sub-functionality 
B may be realized without the use of a Fourier transform, such as by propagation of the 
differently phase changed transformed wavefronts over an extended space. Finally, 
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functionality B requires detection of the intensity characteristics of plurality of 
differently phase changed transformed wavefronts. The outputs of such detection are the 
intensity maps, examples of which are designated by reference numerals 130, 132 and 
134. 

As seen in Fig. 1 A, the third sub-functionality, designated "C" may be realized 
by the following functionalities: 

expressing, such as by employing a computer 136, the plurality of intensity 
maps, such as maps 130, 132 and 134, as at least one mathematical function of phase 
and amplitude of the wavefront being analyzed and of the plurality of different phase 
changes, wherein at least one and possibly both of the phase and the amplitude are 
unknown and the plurality of different phase changes, typically represented by optical 
path delays 1 10, 1 12 and 1 14 to the transformed wavefront 108, are known; and 

employing, such as by means of the computer 136, the at least one 
mathematical function to obtain an indication of at least one and possibly both of the 
phase and the amplitude of the wavefront being analyzed, here represented by the phase 
function designated by reference numeral 138 and the cimplitude function designated by 
reference numeral 139, which, as can be seen, re^spectively represent the phase 
characteristics 102 and the amplitude characteristics 103 of the wavefront 100. In this 
example, wavefront 1 00 may represent the information contained in the compact disk or 
DVD 104, 

In accordance with an embodiment of the present invention, the plurality of 
intensity maps comprises at least four intensity maps. In such a case, employing the 
plurality of intensity maps to obtain an output indicating at least the phase of the 
wavefront being analyzed includes employing a plurality of combinations, each of at 
least three of the plurality of intensity maps, to provide a plurality of indications at least 
of the phase of the wavefront being analyzed. 

Preferably, the methodology also include;s employing the plurality of 
indications of at least the phase of the wavefront being analyzed to provide an enhanced 
indication at least of the phase of the wavefront being analyzed. 

Also in accordance with an embodiment of tlie present invention, the plurality 
of intensity maps comprises at least four intensity majDS. In such a case, employing the 
plurality of intensity maps to obtain an output indicating at least the amplitude of the 
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wavefront being analyzed includes employing a plurality of combinations, each of at 
least three of the plurality of intensity maps, to provide a plurality of indications at least 
of the amplitude of the wavefront being analyzed. 

Preferably, the methodology also include!? employing the plurality of 
indications of at least the amplitude of the wavefront being analyzed to provide an 
enhanced indication at least of the amplitude of the wavefront being analyzed. 

It is appreciated that in this manner, enhanced indications of both phase and 
amplitude of the wavefront may be obtained. 

In accordance with a preferred embodiment of the present invention, at least 
some of the plurality of indications of the amplitude and phase are at least second order 
indications of the amplitude and phase of the wavefront being analyzed. 

In accordance with one preferred embodiment of the present invention, the 
plurality of intensity maps are employed to provide an analytical output indicating the 
amplitude and phase. 

Preferably, the phase changed transformed wavefronts are obtained by 
interference of the wavefront being analyzed along a common optical path. 

In accordance with one preferred embodiment of the present invention, the 
plurality of differently phase changed transformed wavefronts are realized in a manner 
substantially different from performing a delta-function phase change to the transformed 
wavefront, whereby a delta-function phase change is ajDplying a uniform phase delay to 
a small spatial region , having the characteristics of a cLelta-ftmction, of the transformed 
wavefront. 

In accordance with another preferred embodiment of the present invention, the 
plurality of intensity maps are employed to obtain an output indicating the phase of the 
wavefront being analyzed, which is substantially free from halo and shading off 
distortions, which are characteristic of many of the existing 'phase-contrast' methods. 

In accordance with another embodiment of the present invention the output 
indicating the phase of the wavefront being analyzed may be processed to obtain the 
polarization mode of the wavefront being analyzed. 

In accordance with still another embodiment of the present invention, the 
plurality of intensity maps may be employed to obtain an output indicating the phase of 
the wavefront being analyzed by combining the plurality of intensity maps into a second 
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plurality of combined intensity maps, the second plurality being less than the first 
plurality, obtaining at least an output indicative of the phase of the wavefront being 
analyzed from each of the second plurality of combined intensity maps and combining 
the outputs to provide an enhanced indication of the phase of the wavefront being 
analyzed. 

In accordance with yet another embodiment of the present invention, the 
plurality of intensity maps may be employed to obtain an output indicating amplitude of 
the wavefront being analyzed by combining the plurality of intensity maps into a second 
plurality of combined intensity maps, the second plurality being less than the first 
plurality, obtaining at least an output indicative of the amplitude of the wavefront being 
analyzed from each of the second plurality of combined intensity maps and combining 
the outputs to provide an enhanced indication of the amplitude of the wavefront being 
analyzed. 

Additionally in accordance with a preferred embodiment of the present 
invention, the foregoing methodology may be employed for obtaining a plurality of 
differently phase changed transformed wavefronts con-esponding to a wavefront being 
analyzed, obtaining a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts and employing the plurality of intensity maps to obtain an 
output of an at least second order indication of phase of the wavefront being analyzed. 

Additionally or alternatively in accordance with a preferred embodiment of the 
present invention, the foregoing methodology may be employed for obtaining a plurality 
of differently phase changed transformed wavefronts; corresponding to a wavefront 
being analyzed, obtaining a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts and employing the plurality of intensity maps to obtain an 
output of an at least second order indication of amplitude of the wavefront being 
analyzed. 

In accordance with yet another embodiment of the present invention, the 
obtaining of the plurality of differently phase ch.anged transformed wavefronts 
comprises applying a transform to the wavefront being analyzed, thereby to obtain a 
transformed wavefront, and then applying a plurality of different phase and amplitude 
changes to the transformed wavefront, where each of these changes can be a phase 
change, an amplitude change or a combined phase and amplitude change, thereby to 
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obtain a plurality of differently phase and amplitude changed transformed wavefronts. 

In accordance with yet another embodiment of the present invention, a 
wavefront being analyzed comprises at least two wavelength components. In such a 
case, obtaining a plurality of intensity maps also includes dividing the phase changed 
transformed wavefronts according to the at least two wavelength components in order to 
obtain at least two wavelength components of the phase changed transformed 
wavefronts and in order to obtain at least two sets of intensity maps, each set 
corresponding to a different one of the at least two wavelength components of the phase 
changed transformed wavefronts. 

Subsequently, the plurality of intensity maps are employed to provide an output 
indicating the amplitude and phase of the wavefront being analyzed by obtaining an 
output indicative of the phase of the wavefront being analyzed from each of the at least 
two sets of intensity maps and combining the outputs to provide an enhanced indication 
of phase of the wavefront being analyzed. In the enhanced indication, there is no 2n 
ambiguity once the value of the phase exceeds 27u, which conventionally results when 
detecting a phase of a single wavelength wavefront. 

It is appreciated that the wavefront being analyzed may be an acoustic radiation 
wavefront. 

It is also appreciated that the wavefront being analyzed may be an 
electromagnetic radiation wavefront, of any suitable w^avelength, such as visible light, 
infrared, ultra-violet and X-ray radiation. 

It is further appreciated that wavefront 100 may be represented by a relatively 
small number of point sources and defined over a relatively small spatial region. In such 
a case, the detection of the intensity characteristics of the plurality of differently phase 
changed transformed wavefronts may be performed by a detector comprising only a 
single detection pixel or several detection pixels. Additionally, the output indicating at 
least one and possibly both of the phase and amplitude of the wavefront being analyzed, 
may be provided by computer 136 in a straight- forward manner. 

Reference is now made to Fig, IB, which is a simplified partially schematic, 
partially block diagram illustration of a wavefront analysis system suitable for carrying 
out the functionality of Fig. lA in accordance with a preferred embodiment of the 
present invention. As seen in Fig. IB, a wavefront, here designated by reference 
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numeral 150 is focused, as by a lens 152, onto a phase manipulator 154, which is 
preferably located at the focal plane of lens 152. The phase manipulator 154 generates 
phase changes, and may be, for example, a spatial light modulator or a series of 
different transparent, spatially non-uniform objects. 

A second lens 156 is arranged so as to image wavefront 150 onto a detector 
158, such as a CCD detector. Preferably the second lens 156 is arranged such that the 
detector 158 lies in its focal plane. The output of detector 158 is preferably supplied to 
data storage and processing circuitry 160, which preferably carries out functionality "C" 
described hereinabove with reference to Fig. 1 A. 

Reference is now made to Fig. 2, which is a simplified functional block 
diagram illustration of the functionality of Fig. 1 A where time-varying phase changes 
are applied to a transformed wavefront. As seen in Fig, 2, and as explained hereinabove 
with reference to Fig. lA, a wavefront 200 is preferably transformed to provide a 
transformed wavefront 208. 

A first phase change, preferably a spatial phase change, is applied to the 
transformed wavefront 208 at a first time Tl, as indicated by reference numeral 210, 
thereby producing a phase changed transformed wavefront 212 at time Tl. This phase 
changed transformed wavefront 212 is detected, as by detector 158 (Fig. IB), producing 
an intensity map, an example of which is designated by reference numeral 214, which 
map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a second phase change, preferably a spatial phase change, is applied 
to the transformed wavefront 208 at a second time T2, as indicated by reference 
numeral 220, thereby producing a phase changed translbrmed wavefront 222 at time T2. 
This phase changed transformed wavefront 222 is detected, as by detector 158 (Fig. 
IB), producing an intensity map, an example of which is designated by reference 
numeral 224, which map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a third phase change, preferably a spatial phase change, is applied 
to the transformed wavefront 208 at a third time T3, as indicated by reference numeral 
230, thereby producing a phase changed transformed wavefront 232 at time T3. This 
phase changed transformed wavefront 232 is detected, as by detector 158 (Fig. IB), 
producing an intensity map, an example of which is designated by reference numeral 
234, which map is stored as by circuitry 160 (Fig. IB). 
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It is appreciated that any suitable number of spatial phase changes may be 
made at successive times and stored for use in accordance with the present invention. 

In accordance with a preferred embodiment of the present invention, at least 
some of the phase changes 210, 220 and 230, are spatial phase changes effected by 
applying a spatial phase change to part of the transformed wavefront 208. 

In accordance with another preferred embodiment of the present invention, at 
least some of the phase changes 210, 220 and 230, are spatial phase changes, effected 
by applying a time-varying spatial phase change to part of the transformed wavefront 
208. 

In accordance with another preferred embodiment of the present invention, at 
least some of the phase changes 210, 220 and 230, are spatial phase changes, effected 
by applying a non time-varying spatial phase change to part of transformed wavefront 
208, producing spatially phase changed transformed wavefronts 212, 222 and 232, 
which subsequently produce spatially varying intensity maps 214, 224 and 234 
respectively. 

Reference is now made to Fig. 3, which is a simplified functional block 
diagram illustration of the functionality of Fig. 1 A where time- varying phase changes 
are applied to a wavefront prior to transforming thereojf. As seen in Fig. 3, a first phase 
change, preferably a spatial phase change, is applied to a wavefront 300 at a first time 
Tl, as indicated by reference numeral 310. Following application of the first phase 
change to wavefront 300, a transform, preferably a Fourier transform, is applied thereto, 
thereby producing a phase changed transformed waveftont 312 at time Tl. This phase 
changed transformed wavefront 312 is detected, as by detector 158 (Fig. IB), producing 
an intensity map, an example of which is designated by reference numeral 314, which 
map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a second phase change, preferably a spatial phase change, is applied 
to wavefront 300 at a second time T2, as indicated by reference numeral 320. Following 
application of the second phase change to wavefront 300, a transform, preferably a 
Fourier transform, is applied thereto, thereby producing a phase changed transformed 
wavefront 322 at time T2. This phase changed transformed wavefront 322 is detected, 
as by detector 158 (Fig. IB), producing an intensity^ map, an example of which is 
designated by reference numeral 324, which map is stored as by circuitry 160 (Fig. IB). 
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Thereafter, a third phase change, preferably a spatial phase change, is applied 
to wavefront 300 at a third time T3, as indicated by reference numeral 330. Following 
appUcation of the third phase change to wavefront 300, a transform, preferably a Fourier 
transform, is applied thereto, thereby producing a phase changed transformed wavefront 
332 at time T3. This phase changed transformed wavefront 332 is detected, as by 
detector 158 (Fig. IB), producing an intensity map, an example of which is designated 
by reference numeral 334, which map is stored as by circuitry 160 (Fig. IB). 

It is appreciated that any suitable number of spatial phase changes may be 
made at successive times and stored for use in accordance with the present invention. 

In accordance with a preferred embodiment of the present invention, at least 
some of the phase changes 310, 320 and 330, are spatial phase changes effected by 
applying a spatial phase change to part of wavefront 300. 

In accordance with another preferred embodiment of the present invention, at 
least some of the phase changes 310, 320 and 330, are spatial phase changes, effected 
by applying a time-varying spatial phase change to part of wavefront 300. 

In accordance with another preferred embodiment of the present invention, at 
least some of the phase changes 310, 320 and 330, are spatial phase changes, effected 
by applying a non time-varying spatial phase cha:nge to part of wavefront 300, 
producing spatially phase changed transformed wave fronts 312, 322 and 332, which 
subsequently produce spatially varying intensity maps 314, 324 and 334 respectively. 

Reference is now made to Fig. 4, which is a simplified functional block 
diagram illustration of the functionality of Fig. 2, specifically in a case where 
time-varying, non-spatially varying, spatial phase changes are applied to a transformed 
wavefront. As seen in Fig. 4, and as explained hereinal^ove with reference to Fig. 1 A, a 
wavefront 400 is preferably transformed to provide a transformed wavefront 408. A 
preferred transform is a Fourier transform. 

A first spatial phase change is applied to the transformed wavefront 408 at a 
first time Tl, as indicated by reference numeral 410. This phase change preferably is 
effected by applying a spatially uniform spatial phase delay D, designated by reference 
'D=Dr, to a given spatial region of the transformed wavefront 408. Thus, at the given 
spatial region of the transformed wavefront, the value of the phase delay at time Tl is 
Dl, while at the remainder of the transformed wav<^front, where no phase delay is 
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applied, the value of the phase delay is D=0. 

The first spatial phase change 410 thereby produces a spatially phase changed 
transformed wavefront 412 at time TL This spatially phase changed transformed 
wavefront 412 is detected, as by detector 158 (Fig. IB), producing a spatially varying 
intensity map, an example of which is designated by reference numeral 414, which map 
is stored as by circuitry 160 (Fig. IB). 

Thereafter, a second spatial phase change is applied to the transformed 
wavefront 408 at a second time T2, as indicated by rej'erence numeral 420. This phase 
change preferably is effected by applying a spatially uniform spatial phase delay D, 
designated by reference 'D==D2% to a given spatial region of the transformed wavefront 
408. Thus, at the given spatial region of the transfonned wavefront, the value of the 
phase delay at time T2 is D2, while at the remainder of the transformed wavefront, 
where no phase delay is applied, the value of the phase delay is D=0. 

The second spatial phase change 420 thereby produces a spatially phase 
changed transformed wavefront 422 at time T2. This spatially phase changed 
transformed wavefront 422 is detected, as by detector 158 (Fig. IB), producing a 
spatially varying intensity map, an example of which is designated by reference numeral 
424, which map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a third spatial phase change is applied to the transformed wavefront 
408 at a third time T3, as indicated by reference numeral 430. This phase change 
preferably is effected by applying a spatially uniform spatial phase delay D, designated 
by reference 'D=D3', to a given spatial region of the transformed wavefront 408. Thus, 
at the given spatial region of the transformed wavefront, the value of the phase delay at 
time T3 is D3, while at the remainder of the transformed wavefront, where no phase 
delay is applied, the value of the phase delay is D=0. 

The third spatial phase change 430 thereby produces a spatially phase changed 
transformed wavefront 432 at time T3. This spatially phase changed transformed 
wavefront 432 is detected, as by detector 158 (Fig. IB), producing a spatially varying 
intensity map, an example of which is designated by reference numeral 434, which map 
is stored as by circuitry 160 (Fig. IB). 

It is appreciated that any suitable number of spatial phase changes may be 
made at successive times and stored for use in accordance with the present invention. 
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In accordance with a preferred embodiment of the present invention, the 
transform applied to the wavefront 400 is a Fourier transform, thereby providing a 
Fourier-transformed wavefront 408. In addition, the plurality of phase changed 
transformed wavefronts 412, 422 and 432 may be further transformed, preferably by a 
Fourier transform, prior to detection thereof. 

In accordance with a preferred embodiment of the present invention, the spatial 
region of the transformed wavefront 408 to which the spatially uniform, spatial phase 
delays Dl, D2 and D3 are applied at times Tl, T2 and T3 respectively is a spatially 
central region of the transformed wavefront 408. 

In accordance with an embodiment of the present invention, a phase 
component comprising relatively high frequency components may be added to the 
wavefront 400 prior to applying the transform thereto, in order to increase the 
high-frequency content of the transformed wavefront 408 prior to applying the spatially 
uniform, spatial phase delays to a spatial region thereof 

Additionally, in accordance with a preferred embodiment of the present 
invention, the spatial region of the transformed wavefront 408 to which the spatially 
uniform, spatial phase delays Dl, D2 and D3 are applied at times Tl, T2 and T3 
respectively is a spatially central region of the transformed wavefront 408, the transform 
applied to the wavefront 400 is a Fourier transform, and the plurality of phase changed 
transformed wavefronts 412, 422 and 432 are Fourier transformed prior to detection 
thereof. 

In accordance with another embodiment of the^ present invention, the region of 
the transformed wavefront 408 to which the spatially uniform, spatial phase delays Dl, 
D2 and D3 are applied at times Tl, T2 and T3 respectively is a spatially centered 
generally circular region of the transformed wavefront 408. 

In accordance with yet another embodiment of the present invention, the region 
of the transformed wavefront 408 to which the spatially uniform, spatial phase delays 
Dl, D2 and D3 are applied at times Tl, T2 and T3 respectively is a region covering 
approximately one half of the entire region in which transformed wavefront 408 is 
defined. 

In accordance with a preferred embodiment of the present invention, the 
transformed wavefront 408 includes a non-spatially modulated region, termed a DC 
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region, which represents an image of a light source generating the wavefront 400, and a 
non-DC region. The region of the transformed wavefront 408 to which the spatially 
uniform, spatial phase delays DL D2 and D3 are applied at times Tl, T2 and T3 
respectively includes at least parts of both the DC region and the non-DC region. 

Reference is now made to Fig, 5, which is a simplified functional block 
diagram illustration of the functionality of Fig. 3, where time-varying, non-spatially 
varying, spatial phase changes are applied to a wavefront prior to transforming thereof. 

As seen in Fig, 5, a first spatial phase change is appHed to a wavefront 500 at a 
first time Tl, as indicated by reference numeral 510. This phase change preferably is 
effected by applying a spatially uniform spatial phase delay D, designated by reference 
'D=Dr, to a given spatial region of the wavefront 500. Thus, at the given spatial region 
of the wavefront, the value of the phase delay at time Tl is Dl, while at the remainder 
of the wavefront, where no phase delay is applied, the value of the phase delay is D=0. 

Following application of the first spatial phase change to wavefi*ont 500, a 
transform, preferably a Fourier transform, is applied thereto, thereby producing a 
spatially phase changed transformed wavefront 512 at time TL This spatially phase 
changed transformed wavefront 512 is detected, as by detector 158 (Fig. IB), producing 
a spatially varying intensity map, an example of which is designated by reference 
numeral 514, which map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a second spatial phase change is^ applied to wavefront 500 at a 
second time T2, as indicated by reference numeral 520. This phase change preferably is 
effected by applying a spatially uniform spatial phase delay D, designated by reference 
'D=D2', to a given spatial region of the wavefront 500. Thus, at the given spatial region 
of the wavefront, the value of the phase delay at time T2 is D2, while at the remainder 
of the wavefront, where no phase delay is applied, the value of the phase delay is D=0, 

Following application of the second spatial phase change to wavefront 500, a 
transform, preferably a Fourier transform, is applied thereto, thereby producing a 
spatially phase changed transformed wavefront 522 at time T2. This spatially phase 
changed transformed wavefront 522 is detected, as by detector 158 (Fig. IB), producing 
a spatially varying intensity map, an example of v/hich is designated by reference 
numeral 524, which map is stored as by circuitry 160 (Fig. IB). 

Thereafter, a third spatial phase change is applied to wavefront 500 at a third 



46 



time T3, as indicated by reference numeral 530. This phase change preferably is 
effected by applying a spatially uniform spatial phase delay D, designated by reference 
'D=D3', to a given spatial region of the wavefront 500. Thus, at the given spatial region 
of the wavefront, the value of the phase delay at time T3 is D3, while at the remainder 
of the wavefront, where no phase delay is applied, the value of the phase delay is D=0. 

Following application of the third spatial phase change to wavefront 500, a 
transform, preferably a Fourier transform, is applied thereto, thereby producing a 
spatially phase changed transformed wavefront 532 at time T3. This spatially phase 
changed transformed wavefront 532 is detected, as by detector 158 (Fig. IB), producing 
a spatially varying intensity map, an example of wliich is designated by reference 
numeral 534, which map is stored as by circuitry 160 (Fig. IB). 

It is appreciated that any suitable number of spatial phase changes may be 
made at successive times and stored for use in accordance with the present invention. 

In accordance with a preferred embodiment of the present invention, the spatial 
region of the wavefront 500 to which the spatially uniform, spatial phase delays Dl, D2 
and D3 are applied at times Tl, T2 and T3 respectively is a spatially central region of 
the wavefront 500. 

In accordance with an embodiment of the present invention, a phase 
component comprising relatively high frequency coinponents may be added to the 
wavefront 500 prior to applying the spatial phase chaiiges thereto, in order to increase 
the high-frequency content of the wavefront 500. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the spatial region of the wavefront 500 to wliich the spatially uniform, spatial 
phase delays Dl, D2 and D3 are applied at times Tl, T2 and T3 respectively is a 
spatially central region of the wavefront 500, the transiforms are Fourier transforms, and 
the plurality of phase changed transformed wavefronts 512, 522 and 532 are Fourier 
transformed prior to detection thereof. 

In accordance with another embodiment of the present invention, the region of 
the wavefront 500 to which the spatially uniform, spatial phase delays Dl, D2 and D3 
are applied at times Tl, T2 and T3 respectively is a spatially centered generally circular 
region of the wavefront 500. 

In accordance with yet another embodiment of the present invention, the region 
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of the wavefront 500 to which the spatially uniform, spatial phase delays Dl, D2 and 
D3 are applied at times Tl, T2 and T3 respectively is a region covering approximately 
one half of the entire region in which wavefront 500 is defined. 

In accordance with a preferred embodiment of the present invention, the 
wavefront 500 includes a non-spatially modulated region, termed a DC region, which 
represents an image of a light source generating the wavefront 500, and a non-DC 
region. The region of the wavefront 500 to which the spatially uniform, spatial phase 
delays Dl, D2 and D3 are applied at times Tl, T2 and T3 respectively includes at least 
parts of both the DC region and the non-DC region. 

Reference is now made to Fig. 6, which is a simplified functional block 
diagram illustration of the functionality of Fig, lA wheire phase changes are applied to a 
plurality of different wavelength components of a transformed wavefront. As seen in 
Fig. 6, a wavefront 600, which comprises a plurality of different wavelength 
components, is preferably transformed to obtain a transformed wavefront 602. The 
transform is preferably a Fourier transform. 

Similarly to wavefront 600, the transformed wavefront 602 also includes a 
plurality of different wavelength components, represented by reference numerals 604, 
606 and 608. It is appreciated that both the wavefront 600 and the transformed 
wavefront 602 can include any suitable number of wavelength components. 

A plurality of phase changes, preferably spatial phase changes, represented by 
reference numerals 610, 612 and 614 are applied to respective wavelength components 
604, 606 and 608 of the transformed wavefiront, tliereby providing a plurality of 
differently phase changed transformed wavefront com]:>onents, represented by reference 
numerals 620, 622 and 624 respectively. 

The phase changed transformed wavefront components 620, 622, and 624 may 
be transformed, preferably by a Fourier transform, and are subsequently detected, as by 
detector 158 (Fig. IB), producing spatially varying intensity maps, examples of which 
are designated by reference numerals 630, 632 and 634 respectively. These intensity 
maps are subsequently stored as by circuitry 160 (Fig. IB). 

In accordance with an embodiment of the jKesent invention, phase changes 
610, 612 and 614 are effected by passing the transformed wavefront 602 through an 
object, at least one of whose thickness and refracti\^e index varies spatially, thereby 
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applying a different spatial phase delay to each of the wavelength components 604, 606 
and 608 of the transformed wavefront. 

In accordance with another embodiment of the present invention, the phase 
changes 610, 612 and 614 are effected by reflecting the transformed wavefront 602 
from a spatially varying surface, thereby applying a different spatial phase delay to each 
of the wavelength components 604, 606 and 608 of the transformed wavefront. 

In accordance with yet another embodiment of the present invention, the phase 
changes 610, 612 and 614 are realized by passing the transformed wavefront 602 
through a plurality of objects, each characterized in that at least one of its thickness and 
refractive index varies spatially. The spatial variance of the thickness or of the refractive 
index of the plurality of objects is selected in a way such that the phase changes 610, 
612 and 614 differ to a selected predetermined extent for at least some of the plurality of 
different wavelength components 604, 606 and 608. 

Alternatively, the spatial variance of the thickness or refractive index of the 
plurality of objects is selected in a way such that the phase changes 610, 612 and 614 
are identical for at least some of the plurality of different wavelength components 604, 
606 and 608. 

Additionally, in accordance with an embodiment of the present invention, the 
phase changes 610, 612 and 614 are time-varying spatial phase changes. In such a case, 
the plurality of phase changed transformed wavefront components 620, 622 and 624 
include a plurality of differently phase changed transformed wavefronts for each 
wavelength component thereof, and the intensity maps 630, 632 and 634 include a 
time-varying intensity map for each such wavelength component. 

In accordance with an embodiment of the present invention, termed a "white 
light" embodiment, all the wavelength components may be detected by a single 
detector, resulting in a time-varying intensity map representing several wavelength 
components. 

In accordance with another embodiment of the present invention, the plurality 
of phase changed transformed wavefront components 620, 622 and 624 are broken 
down into separate wavelength components, such as by a spatial separation effected, for 
example, by passing the phase changed transformed wavefront components through a 
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dispersion element. In such a case, the intensity maps 630, 632 and 634 are provided 
simultaneously for all of the plurality of different wavelength components. 

Reference is now made to Fig. 7, which is a simplified functional block 
diagram illustration of the functionality of Fig. lA where phase changes are applied to a 
plurality of different wavelength components of a wavefront, prior to transforming 
thereof. As seen in Fig. 7, a wavefront 700 comprises a plurality of different wavelength 
components 704, 706 and 708. It is appreciated that the wavefront can include any 
suitable number of wavelength components. 

A plurality of phase changes, preferably spatial phase changes, represented by 
reference numerals 710, 712 and 714, are applied to the respective wavelength 
components 704, 706 and 708 of the wavefiront. 

Following application of the spatial phase changes to wavefront components 
704, 706 and 708, a transform, preferably a Fourier transform, is applied thereto, 
thereby providing a plurality of different phase changed transformed wavefront 
components, represented by reference numerals 720, 722 and 724 respectively. 

These phase changed transformed wavefront components 720, 722 and 724 are 
subsequently detected, as by detector 158 (Fig. IB), producing spatially varying 
intensity maps, examples of which are designated by reference numerals 730, 732 and 
734, These intensity maps are subsequently stored as by circuitry 160 (Fig. IB). 

In accordance with an embodiment of the present invention, phase changes 
710, 712 and 714 are effected by passing the wavefront 700 through an object, at least 
one of whose thickness and refractive index varies spatially, thereby applying a 
different spatial phase delay to each of the wavelength components 704, 706 and 708 of 
the wavefront. 

In accordance with another embodiment of the present invention, the phase 
changes 710, 712 and 714 are effected by reflecting tlie wavefront 700 from a spatially 
varying surface, thereby applying a different spatial phase delay to each of the 
wavelength components 704, 706 and 708 of the wavefront. 

In accordance with yet another embodiment of the present invention phase 
changes 710, 712 and 714 are realized by passing the wavefront 700 through a plurality 
of objects, each characterized in that at least one of its thickness and refractive index 
varies spatially. The spatial variance of the thickness or refractive index of these objects 



50 



is selected in a way such that the phase changes 710, 712 and 714 differ to a selected 
predetermined extent for at least some of the plu]rality of different wavelength 
components 704, 706 and 708. 

Alternatively, the spatial variance of the thickness or refractive index of these 
objects is selected in a way that the phase changes 710, 712 and 714 are identical for at 
least some of the plurality of different wavelength components 704, 706 and 708. 

Reference is now made to Fig. 8, which is a simplified functional block 
diagram illustration of the functionality of Fig. lA where phase changes are applied to a 
plurality of different polarization components of a transformed wavefront. As seen in 
Fig. 8, a wavefront 800, which comprises a plui'ality of different polarization 
components, is preferably transformed to obtain a transformed wavefront 802. The 
transform is preferably a Fourier transform. Similarly to wavefront 800, the transformed 
wavefront 802 also includes a plurality of different polarization components, 
represented by reference numerals 804 and 806. It is appreciated that the polarization 
components 804 and 806 can be either spatially different or spatially identical, but are 
each of different polarization. It is further appreciated that both the wavefront 800 and 
the transformed wavefront 802 preferably each include two polarization components but 
can include any suitable number of polarization components. 

A plurality of phase changes, preferably spatial phase changes, represented by 
reference numerals 810 and 812, are applied to the respective polarization components 
804 and 806 of the transformed wavefront 802, thereby providing a plurality of 
differently phase changed transformed wavefront components, represented by reference 
numerals 820 and 822 respectively. 

It is appreciated that phase changes 810 and 812 can be different for at least 
some of the plurality of different polarization components 804 and 806. Alternatively, 
phase changes 810 and 812 can be identical for at least some of the plurality of different 
polarization components 804 and 806. 

The phase changed transformed wavefront components 820 and 822 are 
detected, as by detector 158 (Fig. IB), producing spatially varying intensity maps, 
examples of which are designated by reference numerals 830 and 832, These intensity 
maps are subsequently stored as by circuitry 160 (Fig. IB). 

Reference is now made to Fig. 9, which is a simplified functional block 
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diagram illustration of the functionality of Fig. lA whers phase changes are applied to a 
plurality of different polarization components of a w^avefront prior to transforming 
thereof. As seen in Fig. 9, a wavefront 900 comjDrises a plurality of different 
polarization components 904 and 906. It is appreciated that the wavefront preferably 
includes two polarization components but can include any suitable number of 
polarization components, 

A plurality of phase changes, preferably spatial phase changes, represented by 
reference numerals 910 and 912, are applied to the respective polarization components 
904 and 906 of the wavefront. 

It is appreciated that phase changes 910 and 912 can be different for at least 
some of the plurality of different polarization components 904 and 906. Alternatively, 
phase changes 910 and 912 can be set to be identical for at least some of the plurality of 
different polarization components 904 and 906. 

Following application of the spatial phase chimges to wavefront components 
904 and 906, a transform, preferably a Fourier transfiDrm, is applied thereto, thereby 
providing a plurality of different phase changed transformed wavefront components, 
designated by reference numerals 920 and 922 respectively. 

Phase changed transformed wavefront components 920 and 922 are 
subsequently detected, as by detector 158 (Fig. IB), producing spatially varying 
intensity maps, examples of which are designated by reference numeral 930 and 932. 
These intensity maps are subsequently stored as by circuitry 160 (Fig. IB). 

Reference is now made to Fig. lOA, which is a simplified functional block 
diagram illustration of the functionality of Fig. lA, where a wavefront being analyzed 
comprises at least one one-dimensional component. In the embodiment of Fig. lOA, a 
one-dimensional Fourier transform is applied to the wa^^efront. Preferably, the transform 
is performed in a dimension perpendicular to a direction of propagation of the 
wavefront being analyzed, thereby to obtain at least one one-dimensional component of 
the transformed wavefront in the dimension perp^endicular to the direction of 
propagation. 

A plurality of different phase changes are applied to each of the at least one 
one-dimensional components, thereby obtaining at least one one-dimensional 
component of the plurality of phase changed transformed wavefronts. 
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A plurality of intensity maps are employed to obtain an output indicating 
amplitude and phase of the at least one one-dimensional component of the wavefront 
being analyzed. 

As seen in Fig. lOA, a plurality of different phase changes are appUed to at 
least one one-dimensional component of a transformed wavefront. In the illustrated 
embodiment, typically five one-dimensional components of a wavefront are shown and 
designated by reference numerals 1001, 1002, 1003, 1004 and 1005. The wavefront is 
transformed, preferably by a Fourier transform. It is thus appreciated that due to 
transform of the wavefront, the five one-dimensional components 1001, 1002, 1003, 
1004 and 1005 are transformed into five corresponding one-dimensional components of 
the transformed wavefront, respectively designated by reference numerals 1006, 1007, 
1008, 1009 and 1010. 

Three phase changes, respectively designated 1011, 1012 & 1013 are each 
applied to the one-dimensional components 1006, 1007, 1008, 1009 and 1010 of 
transformed wavefront to produce three phase chiuiged transformed wavefronts, 
designated generally by reference numerals 1016, 1018 and 1020. 

In the illustrated embodiment, phase changed transformed wavefi-ont 1016 
includes five one-dimensional components, respectively designated by reference 
numerals 1021, 1022, 1023, 1024 and 1025. 

In the illustrated embodiment, phase changed transformed wavefront 1018 
includes five one-dimensional components, respectively designated by reference 
numerals 1031, 1032, 1033, 1034 and 1035. 

In the illustrated embodiment, phase changed transformed wavefront 1020 
includes five one-dimensional components, respectively designated by reference 
numerals 1041, 1042, 1043, 1044 and 1045. 

The phase changed transformed wavefronts 1016, 1018 and 1020 are detected, 
as by detector 158 (Fig. IB), producing three intensity maps, designated generally by 
reference numerals 1046, 1048 and 1050. 

In the illustrated embodiment, intensity map 1046 includes five 
one-dimensional intensity map components, respectively designated by reference 
numerals 1051, 1052, 1053, 1054 and 1055. 

In the illustrated embodiment, intensity map 1048 includes five 
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one-dimensional intensity map components, respectively designated by reference 
numerals 1061, 1062, 1063, 1064 and 1065. 

In the illustrated embodiment, intensity map 1050 includes five 
one-dimensional intensity map components, respectively designated by reference 
numerals 1071, 1072, 1073, 1074 and 1075. 

The intensity maps 1046, 1048 and 1050 are stored as by circuitry 160 (Fig. 

IB). 

In accordance with an embodiment of the present invention, the wavefront 
being analyzed, illustrated in Fig. lOA by the one-dimensional components 1001, 1002, 
1003, 1004 and 1005, may comprise a plurality of different wavelength components and 
the plurality of different phase changes, 1011, 1012 and 1013, are applied to the 
plurality of different wavelength components of each of the plurality of 
one-dimensional components of the wavefront being a.nalyzed. Preferably, obtaining a 
plurality of intensity maps 1046, 1048 and 1050, includes dividing the plurality of 
one-dimensional components of the plurality of phase changed transformed wavefronts 
1016, 1018 and 1020 into separate wavelength components. 

Preferably, dividing the plxirality of one-dimensional components of the 
plurality of phase changed transformed wavefronts into separate wavelength 
components is achieved by passing the plurality of phase changed transformed 
wavefronts 1016, 1018 and 1020 through a dispersion element. 

Reference is now made to Fig. lOB, which is a simplified partially schematic, 
partially pictorial illustration of a wavefront analysis s}^stem suitable for carrying out the 
functionality of Fig, lOA in accordance with a preferred embodiment of the present 
invention. 

As seen in Fig. lOB, a wavefront, here designated by reference numeral 1080, 
and here including five one-dimensional components 1081, 1082, 1083, 1084 and 1085 
is focused, as by a cyhndrical lens 1086 onto a single axis displaceable phase 
manipulator 1087, which is preferably located at the focal plane of lens 1086. Lens 1086 
preferably produces a one-dimensional Fourier transform of each of the 
one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 along the 
Y-axis, 

As seen in Fig. lOB, the phase manipulator 1087 preferably comprises a 
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multiple local phase delay element, such as a spatially non-uniform transparent object, 
typically including five different phase delay regions, each arranged to apply a phase 
delay to one of the one-dimensional components at a given position of the object along 
an axis, here designated as the X-axis, extending perjDendicularly to the direction of 
propagation of the wavefront along a Z-axis and perpendicular to the axis of the 
transform produced by lens 1086, here designated as the Y-axis. 

A second lens 1088, preferably a cylindrical lens, is arranged so as to image the 
one-dimensional components 1081, 1082, 1083, 1084 and 1085 onto a detector 1089, 
such as a CCD detector. Preferably the second lens 1088 is arranged such that the 
detector 1089 lies in its focal plane. The output of detector 1089 is preferably supplied 
to data storage and processing circuitry 1090, which preferably carries out functionality 
"C" described hereinabove with reference to Fig. 1 A. 

There is provided relative movement between the optical system comprising 
phase manipulator 1087, lenses 1086 and 1088 and detector 1089 and the 
one-dimensional wavefront components 1081, 1082, 1083, 1084 and 1085 along the 
X-axis. This relative movement sequentially matches different phase delay regions with 
different wavefront components, such that preferably each wavefront component passes 
through each phase delay region of the phase manipulator 1087. 

It is a particular feature of the embodiment of Figs. lOA and lOB, that each of 
the one dimensional components of the wavefront is separately processed. Thus, in the 
context of Fig. lOB, it can be seen that the five one-dimensional wavefront components 

1081, 1082, 1083, 1084 and 1085 are each focused by a separate portion of the 
cylindrical lens 1086, are each imaged by a corresponding separate portion of the 
cylindrical lens 1088 and each pass through a distinct region of the phase manipulator 
1087. The images of each of the five one-dimensional wavefront components 1081, 

1082, 1083, 1084 and 1085 at detector 1089 are thus seen to be separate and distinct 
images, as designated respectively by reference numerals 1091, 1092, 1093, 1094 and 
1095. It is appreciated that these images may appeju: on separate detectors together 
constituting detector 1089 instead of on a monolithic detector. 

In accordance with an embodiment of the present invention, the transform 
applied to the wavefront includes an additional Fc^urier transform. This additional 
Fourier transform may be performed by lens 1086 or by an additional lens and is 
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operative to minimize cross-talk between different one-dimensional components of the 
wavefront. In such a case, preferably a further transform is applied to the phase changed 
transformed wavefront. This further transform may be performed by lens 1088 or by an 
additional lens. 

Reference is now made to Fig. 11, which is a simplified functional block 
diagram illustration of the functionality of Fig. 1 A, where an additional transform is 
applied following the application of spatial phase changes. As seen in Fig. 11, and as 
explained hereinabove with reference to Fig, lA, a wavefront 1100 is transformed, 
preferably by a Fourier transform and a plurality of phase changes are applied to the 
transformed wavefront, thereby to provide a plurality of differently phased changed 
transformed wavefronts, represented by reference numerals 1 120, 1 122, and 1 124. 

The phase changed transformed wavefronts are subsequently transformed, 
preferably by a Fourier transform, and then detected, as by detector 158 (Fig. IB), 
producing spatially varying intensity maps, examples of which are designated by 
reference numerals 1130, 1132 and 1134. These intensity maps are subsequently stored 
as by circuitry 160 (Fig. IB). 

It is appreciated that any suitable number of differently phased changed 
transformed wavefronts can be obtained, and subsequently transformed to a 
corresponding plurality of intensity maps to be storecl for use in accordance with the 
present invention. 

Reference in now made to Fig. 12, which Is a simplified functional block 
diagram iUustration of the functionality of Fig. lA, wherein intensity maps are 
employed to provide information about a wavefront being analyzed, such as indications 
of amplitude and phase of the wavefront. As seen in Fig. 12, and as explained 
hereinabove with reference to Fig. 1 A, a wavefront 1200 is transformed, preferably by a 
Fourier transform, and phase changed by a phase-ch.ange function to obtain several, 
preferably at least three, differently phase-changed transformed wavefronts, respectively 
designated by reference numerals 1210, 1212 and 1214. The phase changed transformed 
wavefronts 1210, 1212 and 1214 are subsequently detected, as by detector 158 (Fig. 
IB), producing spatially varying intensity maps, examples of which are designated by 
reference numerals 1220, 1222 and 1224. 

In parallel to producing the plurality of intensity maps, such as intensity maps 
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1220, 1222 and 1224, the expected intensity maps are expressed as a first function of 
the amplitude of wavefront 1200, of the phase of wavefront 1200, and of the phase 
change function characterizing the differently phase changed transformed wavefronts 
1210, 1212 and 1214, as indicated at reference numeral 1230. 

In accordance with a preferred embodiment of the present invention, at least 
one of the phase and the amplitude of the wavefront is unknown or both the phase and 
the amplitude are unknown. The phase-change function is known. 

The first function of the phase and amplitude of the wavefront and of the phase 
change function is subsequently solved as indicated at reference numeral 1235, such as 
by means of a computer 136 (Fig. 1 A), resulting in an expression of at least one and 
possibly both of the amplitude and phase of wavefront 1200 as a second function of the 
intensity maps 1220, 1222 and 1224, as indicated at reference numeral 1240. 

The second function is then processed together with the intensity maps 1220, 
1222 and 1224 as indicated at reference numeral 1242. As part of this processing, 
detected intensity maps 1220, 1222 and 1224 are substituted into the second function. 
The processing may be carried out by means of a computer 136 (Fig. lA) and provides 
information regarding wavefi-ont 1200, such as indications of at least one and possibly 
both of the amplitude and the phase of the wavefront. 

In accordance with a further embodiment of th.e present invention, the plurality 
of intensity maps comprises at least four intensity maps. In such a case, employing the 
plurality of intensity maps to obtain an indication of at least one of the phase and the 
amplitude of the wavefront 1200 includes employing a plurality of combinations, each 
of the combinations being a combination of at least ttiree of the plurality of intensity 
maps, to provide a plurality of indications of at least one of the phase and the amplitude 
of wavefront 1200. Preferably, this methodology also includes employing the plurality 
of indications of at least one of the phase and the amplitude of the wavefront 1200 to 
provide an enhanced indication at least one of the phase and the amplitude of the 
wavefront 1200. 

In accordance with a preferred embodiment of the present invention, at least 
some of the plurality of indications of the amplitude and phase are at least second order 
indications of the amplitude and phase of the wavefront 1200. 

In accordance with another embodiment of the present invention, the first 
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function may be solved as a function of some unknowis to obtain the second function 
by expressing, as indicated by reference numeral 1240, some unknowns, such as at least 
one of the amplitude and phase of wavefront 1200, as a second function of the intensity 
maps. 

Accordingly, solving the first function may include: 

defining a complex function of the amplitude of wavefront 1200, of 
the phase of wavefront 1200, and of the phase change function characterizing the 
differently phase changed transformed wavefronts 1210, 1212 and 1214. This complex 
function is characterized in that intensity at each location in the plurality of intensity 
maps is a function predominantly of a value of the complex function at that location and 
of the amplitude and the phase of wavefront 1200 at the same location; 

expressing the complex function as a third function of the plurality of 
intensity maps 1220, 1222 and 1224; and 

obtaining values for the unknowns, such as at least one of phase and 
amplitude of wavefront 1200, by employing the complex function expressed as a 
function of the plurality of intensity maps. 

In accordance with this embodiment, preferably the complex function is a 
convolution of another complex function, which has an amplitude and phase identical to 
the amplitude and phase of wavefront 1200, and of a Fourier transform of the phase 
change function characterizing the differently phase changed transformed wavefronts 
1210, 1212 and 1214. 

Reference in now made to Fig. 13, which is a simplified functional block 
diagram illustration of part of the functionality of Fig. lA, wherein the transform 
applied to the wavefront being analyzed is a Fourier transform, wherein at least three 
different spatial phase changes are applied to the thus transformed wavefront, and 
wherein at least three intensity maps are employed to obtain indications of at least one 
of the phase and the amplitude of the wavefront. 

As explained hereinabove with reference to Fig. lA, a wavefront 100 (Fig. lA) 
being analyzed, is transformed and phase changed by at least three different spatial 
phase changes, all governed by a spatial function, to obtain at least three differently 
phase-changed transformed wavefronts, represented by reference numerals 120, 122 and 
124 (Fig, lA) which are subsequently detected, as by detector 158 (Fig. IB), producing 
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spatially varying intensity maps, examples of which are designated by reference 
numerals 130, 132 and 134 (Fig. lA). As seen in Fig. 13, and designated as 
sub-functionality "C" hereinabove with reference in Fig. lA, the intensity maps are 
employed to obtain an output indication of at least one and possibly both of the phase 
and the amplitude of the wavefront being analyzed. 

Turning to Fig. 13, it is seen that the wavefront being analyzed is expressed as 

a first complex function f(x) = A{x)q'^^''\ where 'x' is a general indication of a spatial 

location. The complex function has an amplitude distribution A(x) and a phase 
distribution (p(x) identical to the amplitude and phase of the wavefront being analyzed. 

The first complex function f(x) = A(x)g'^^''^ is indicated by reference numeral 1300. 

As noted hereinabove with reference to Fig. lA, each of the plurality of 
different spatial phase changes is applied to the transformed wavefront preferably by 
applying a spatially uniform spatial phase delay having a known value to a given spatial 
region of the transformed wavefront. As seen in Fig. 13, the spatial function governing 
these different phase changes is designated by 'G' and an example of which, for a phase 
delay value of 0, is designated by reference numeral 1304. 

Function 'G' is a spatial function of the phase change applied in each spatial 
location of the transformed wavefront. In the specific example designated by reference 
numeral 1304, the spatially uniform spatial phase delay, having a value of 0, is applied 
to a spatially central region of the transformed wavefront, as indicated by the central 
part of the function having a value of 0, which is greater than the value of the function 
elsewhere. 

A plurality of expected intensity maps, indicated by spatial functions Ii(x), l2(x) 
and /i(xj, are each expressed as a function of the first complex function f(x) and of the 
spatial function G, as indicated by reference numeral 1308. 

Subsequently, a second complex function S(x)^ which has an absolute value 
\S(x) I and a phase a(x), is defined as a convolution of the first complex function f(x) and 
of a Fourier transform of the spatial function 'G\ This second complex function, 
designated by reference numeral 1312, is indicated by the equation 

5(x)= /(x)^S(G) = where the symbol indicates convolution and S(g) 

is the Fourier transform of the function 'G'. 

59 



The difference between q)(x), the phase of the wavefront, and a(x), the phase of 
the second complex function, is indicated by y/(x), as designated by reference numeral 
1316. 

The expression of each of the expected intensity maps as a function off(x) and 
G, as indicated by reference numeral 1308, the definition of the absolute value and the 
phase of S(x), as indicated by reference numeral 1312 and the definition of y/(x), as 
indicated by reference numeral 1316, enables expression of each of the expected 
intensity maps as a third function of the amplitude of the wavefront A(x), the absolute 
value of the second complex function \S(x)\, the difference between the phase of the 
wavefront and the phase of the second complex function ii/(x), and the knovm phase 
delay produced by one of the at least three different phase changes which each 
correspond to one of the at least three intensity maps. 

This third function is designated by reference numeral 1320 and includes three 
functions, each preferably having the general form 

2 
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where In(x) are the expected intensity maps and 



n = 1,2 or 3, In the three functions, 6j, 62 and O3 are the known values of the uniform 
spatial phase delays, each applied to a spatial region of the transformed wavefront, thus 
effecting the plurality of different spatial phase changes which produce the intensity 
maps Ii(x), l2(x) and l3(x), respectively. 

It is appreciated that preferably the third fimction at any given spatial location 
xq is a function of A, ^and \S\ only at the same spatial location xq. 

The intensity maps are designated by reference numeral 1324. 

The third function is solved for each of the specific spatial locations xq, by 
solving at least three equations, relating to at least three intensity values Ii(xo)^ hi^) and 
h(xo) at at least three different phase delays 6u Oi and Gs, thereby to obtain at least part 
of three unknowns A(xq), \S(xq)\ and \j/(xo). This process is typically repeated for all 
spatial locations and results in obtaining the amplitude of the wavefi-ont A(x), the 
absolute value of the second complex function \S(x)\ and the difference between the 
phase of the wavefront and the phase of the second complex function xj/(x), as indicated 
by reference numeral 1328. 
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Thereafter, once A(x), \S(x)\ and i//(x) are known, the equation defining the 
second complex function, represented by reference numeral 1312, is typically solved 
globally for a substantial number of spatial locations 'x' to obtain a(x), the phase of the 
second complex function, as designated by reference numeral 1332. 

Finally, the phase (p(x) of the wavefront being analyzed is obtained by adding 
the phase a(x) of the second complex function to the difference y/(x) between the phase 
of the wavefront and the phase of the second complex function, as indicated by 
reference numeral 1336. 

In accordance with an embodiment of the present invention, the absolute value 
151 of the second complex function is obtained preferably for every specific spatial 
location xo by approximating the absolute value to a polynomial of a given degree in the 
spatial location x. 

In accordance with another preferred embodiment of the present invention, the 
phase a(x) of the second complex function is obtained by expressing the second 
complex function S(x) as an eigen-value problem, such as S = S M where M is a 
matrix, and the complex function is an eigen-vector of the matrix obtained by an 
iterative process. An example of such an iterative process is 
5,=|5|, S„,,-5„M/||5„M||, where n is the iterative step number. 

In accordance with yet another preferred embodiment of the present invention, 
the phase a(x) of the second complex function is obtained by approximating the Fourier 
transform of the spatial function 'G% governing the spatial phase change, to a 
polynomial in the location x, by approximating the second complex function S(x) to a 
polynomial in the location x, and by solving, according to these approximations, the 



equation defining the second complex function: S{x) = — — S{x) 

I ^(^) I 
V / 



* S[G] , where 



A(x)e 



the function — VVri — is known. 

In accordance with still another preferred embodiment of the present invention, 
at any location x the amplitude A(x) of the wavefront being analyzed, the absolute value 
\S{x)\ of the second complex function, and the difference i//(x) between the phase of the 
second complex function and the phase of the wavefront are obtained by a best-fit 
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method, such as a least-square method, preferably a linear least-square method, from the 
values of the intensity maps at this location In(x), where n=},2,..,,N and N is the number 
of intensity maps. The accuracy of this process increases as the number N of the 
plurality of intensity maps increases. 

In accordance with one preferred embodiment of the present invention, the 
plurality of different phase changes comprises at least four different phase changes, the 
plurality of intensity maps comprises at least four intensity maps, and the function 
designated by reference numeral 1320 can express each of the expected intensity maps 
as a third function of: 

the amplitude of the wavefront^fxj; 

the absolute value of the second complex function \S(x)\; 

the difference between the phase of the wavefront and the phase of the 
second complex function y/(x); 

the known phase delay produced by one of the at least four different 
phase changes each of which corresponds to one of the at least four intensity maps; and 

at least one additional unknovm relating to the wavefront analysis, 
where the number of the at least one additional unknown is no greater than the number 
by which the plurality intensity maps exceeds three. 

The third function 1320, is then solved by solving at least four equations, 
resulting from at least four intensity values at at least four different phase delays, 
thereby to obtain the amplitude of the wavefront being analyzed, the absolute value of 
the second complex function, the difference between the phase of the wavefront and the 
phase of the second complex function and the at least one additional unknown. 

In accordance with another preferred embodiment of the present invention, the 
values of the uniform spatial phase delays 0^, 02, 6iv applied to a spatial region of the 
transformed wavefront, thus effecting the plurality of different spatial phase changes, 
producing the intensity maps Ii(x), l2(x), In(x) respectively, are chosen as to 
maximize contrast in the intensity maps and to minimize effects of noise on the phase of 
the wavefront being analyzed. 

In accordance with one more preferred embodiment of the present invention, 
the function designated by reference numeral 1320, expressing each of the expected 
intensity maps as a third function of the amplitude of the wavefront A(x), the absolute 
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value of the second complex function \S(x)l the difference between the phase of the 
wavefront and the phase of the second complex function i//(x), and the known phase 
delay 6, produced by one of the at least three different phase changes which each 
correspond to one of the at least three intensity maps, comprises several fiinctionalities: 

defining fourth, fifth and sixth complex functions, designated as /3o(x), f3s(x) 
and /3c(x) respectively, none of which is a function of any of the plurality of intensity 
maps or of the spatial function 'G' governing the phase change. Each of the fourth, fifth 
and sixth complex functions is preferably a function of the amplitude of the wavefront 
A(x). the absolute value of the second complex function \S(x)i the difference between 
the phase of the wavefront and the phase of the second complex function y/(x)\ and 

expressing each of the plurality of intensity maps In(x) as 
/^^(x) = p^{x) + ^^(x)cos(^^) + /?5(x)sin(^„) , where On is the value of the phase delay 
corresponding to intensity map In(x). Each intensity map In(x), where n = 1,2, ... A^, 



preferably expressed as / {x) 



A(x) + (e'^^ -l)\S{x)\e '^^""^ 



can be subsequently 



expressed as (x) = J3q (x) + /J^ W ^^K^n ) + A (^) sin(^n ) ^ where 



'fio(x)=A{xf+2\S(x)\^ -2A(x)\S(x)\cos(y/) 

. (x) = 2A(x) \S(x)\ cos(^) - 2 1 p 
J3,(x) ^2A{x)\Six)\sin{i//) 

Preferably the foregoing methodology also includes solving the third function 
1320 by using a linear least-square method to compute from the different intensities 
the values of j3o, J3c and A best fitting to -P^^P^ + sin 6^ . 

Subsequently the amplitude A(x) is found by A(x) = V>^o W + W ^ ^1^^ absolute 
value \S(x)\ of the second complex function is found by solving the second degree 

2 2 

equation \S{x)t -l3^{x)\S{x)f ^ ^^^""^ ^^^'^""^ =0 for \S(x)t and xKx) is 

2 

found by Vi^^ = argC/^c (^) + 2 1 ^(x) | (^)) 
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In accordance with yet another preferred embodiment of the present invention, 
solving of the third function, designated by reference numeral 1320, to obtain, as 
designated by reference numeral 1328, the amplitude of the wavefront ^(3cj, the absolute 
value of the second complex function \S(x)\ and the difference between the phase of the 
wavefront and the phase of the second complex function y/(x), includes several 
functionalities: 

obtaining two solutions for the absolute value \S(x)\ of the second complex 
function, these two solutions, being designated by \Sh(x)\ and \Si(x)l namely a higher 
value solution and a lower value solution respectively; and 

combining the two solutions into an enhanced absolute value solution \S(x)\ 
for the absolute value of the second complex function, by choosing at each spatial 
location 'xq' either the higher value solution \Sh(xo)\ or the lower value solution \Si(xo)\ 
such that the enhanced absolute value solution satisfies the second complex function, 
designated by reference numeral 1312. 

Preferably the methodology also includes: 

obtaining two solutions for each of the amplitude A(x) of the wavefront 
being analyzed and the difference i//(x) between the phase of the wavefront and the 
phase of the second complex function, these two solutions being higher value solutions 
Ah(x) and y/h(x) and lower value solutions ^/fxj and i//i(x); and 

combining the two solutions Ah(x) and Ai(x) for the amplitude into an 
enhanced amplitude solution A(x) by choosing at each spatial location 'xo' either the 
higher value solution Ah(xo) or the lower value solution Ai(xo) in a way that at each 
spatial location 'xq' if \Sh(xo)\ is chosen for the absolute value solution, then Ah(xQ) is 
chosen for the amplitude solution and at each location 'xf if \Si(xi)\ is chosen for the 
absolute value solution, ih^n Ai(xi) is chosen for the amplitude solution; and 

combining the two solutions y/h(x) and i//i(x) of the difference between 
the phase of the wavefront and the phase of the second complex function into an 
enhanced difference solution y/(x), by choosing at each spatial location 'xq' either the 
higher value solution ii/h(xo) or the lower value solution \//i(xo) in a way that at each 
spatial location 'xq' if \Sh(xo)\ is chosen for the absolute value solution, then y/hfxo) is 
chosen for the difference solution and at each location 'xj' if \Si(xj)\ is chosen for the 
absolute value solution, then y/i(x/) is chosen for the difference solution. 
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Additionally, in accordance with an embodiment of the present invention, the 
plurality of different phase changes applied to the transformed wavefront, thereby to 
obtain a plurality of differently phase changed transformed wavefronts, also include 
amplitude changes, resulting in a plurality of differently phase and amplitude changed 
transformed wavefronts. These amplitude changes are preferably known amplitude 
attenuations applied to the same spatial region of the transformed wavefront to which 
the uniform phase delays 6i, 9^.,,, On, are applied, the spatial region being defined by 
the spatial function 'G\ 

The amplitude attenuations are designated by au cr2,..., a^^ where the n-th 
change, where n = 1,2, ... N, applied to the transformed wavefront includes a phase 
change On and an amplitude attenuation cr« . It is appreciated that some of the phase 
changes may be equal to zero, indicating no phase-change and that some of the 
amplitude attenuations may be equal to unity, indicating no amplitude attenuation. 

In this embodiment, the function designated by reference numeral 1320, 
expressing each of the expected intensity maps In(x) as a third function of the amplitude 
of the wavefront A(x), the absolute value of the second complex function \S(x)\, the 
difference between the phase of the wavefront and the phase of the second complex 
function y/(x), and the phase delay On , also expresses each of the expected intensity 
maps also as a function of the amplitude attenuation cxn and comprises several 
functionalities: 

defining fourth, fifth, sixth and seventh complex functions, designated by 
Po(x), J3j(x), f32(x) and {3s(x) respectively, none of which is a function of any of the 
plurality of intensity maps or of the spatial function 'G' governing the phase and 
amplitude changes. Each of the fourth, fifth, sixth and seventh complex functions is 
preferably a function of the amplitude of the wavefront A(x), the absolute value of the 
second complex function \S(x)l the difference between the phase of the wavefront and 
the phase of the second complex function y/(x)\ 

defining an eighth function, designated as a combination of the phase 
delay and of the amplitude attenuation, where for the n-Xh change applied to the 
transformed wavefront, including a phase change On and an amplitude attenuation an , 
this eighth function is designated by Preferably the combination jUn is defined by 
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expressing each of the pluraUty of intensity maps In(x) as 

fi,{x)^A{x)\S(x)\e-^'^' and /S,(x) ^ A(x)\S(x)\e^^'^^ . 

Preferably the foregoing methodology also includes solving the third function 
1320 by computing from the different intensities In(x), the values of /3o(x), Pi(x), p2(x) 

and f33(x) best fitting to the equation - y9o(x) + + + A W/^n • 

Subsequently the amplitude A(x) is found by A{x) = ^f3^{x) , the absolute value \S(x)\ 

of the second complex function is found by 1 S{x) |- ™Ld y/(x) is found by 

solving e^^*^^^ = cingleip ^{x)) . 

It is appreciated that the amplitude attenuations cu 0-2, > o^n, may be 
unknown. In such a case, additional intensity maps are obtained, where the number of 
the unknowns is no greater than the number by which the plurality of intensity maps 
exceeds three. The unknowns are obtained in a manner similar to that described 
hereinabove, where there exists at least one unknown relating to the wavefront analysis. 

Reference is now made to Fig, 14, which is a simplified partially schematic, 
partially pictorial illustration of part of one preferred embodiment of a wavefront 
analysis system of the type shown in Fig, IB. As seen in Fig. 14, a wavefront, here 
designated by reference numeral 1400 is partially transmitted through a beam splitter 
1402 and subsequently focused, as by a lens 1404 onto a phase manipulator 1406, which 
is preferably located at the focal plane of lens 1404. The phase manipulator 1406 may 
be, for example, a spatial light modulator or a series of different transparent, spatially 
non-uniform objects. 

A reflecting surface 1408 is arranged so as to reflect wavefront 1400 after it 
passes through the phase manipulator 1406. The reflected wavefront is imaged by lens 
1404 onto a detector 1410, such as a CCD detector via beam splitter 1402. Preferably 
the beam splitter 1402 and the detector 1410 are arranged such that the detector 1410 
lies in the focal plane of lens 1404. The output of detector 1410 is preferably supplied to 
data storage and processing circuitry 1412, which preferably carries out functionality 
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"C" described hereinabove with reference to Fig. 1 A. 

It is appreciated that adding the reflecting surface 1408 to an imaging system, 
doubles the phase delay generated by phase manipulator 1406, enables imaging with a 
single lens 1404, and generally enables realization of a more compact system. 

Reference is now made to Fig. 15, which is a simplified partially schematic, 
partially pictorial illustration of a system for surface mapping employing the 
functionality and structure of Figs. lA and IB. As seen in Fig. 15, a beam of radiation, 
such as light or acoustic energy, is supplied from a radiation source 1500, optionally via 
a beam expander 1502, onto a beam splitter 1504, which reflects at least part of the 
radiation onto a surface 1506 to be inspected. The radiation reflected from the inspected 
surface 1506, is a surface mapping wavefiront, which has an amplitude and a phase, and 
which contains information about the surface 1506. At least part of the radiation 
incident on surface 1506 is reflected from the surface 1506 and transmitted via the beam 
splitter 1504 and focused via a focusing lens 1508 onto a phase manipulator 1510, 
which is preferably located at the image plane of radiation source 1500. 

The phase manipulator 1510 may be, for example, a spatial light modulator or 
a series of different transparent, spatially non-uniform objects. It is appreciated that 
phase manipulator 1510 can be configured such that a substantial part of the radiation 
focused thereonto is reflected therefrom. Alternatively the phase manipulator 1510 can 
be configured such that a substantial part of the radiation focused thereonto is 
transmitted therethrough. 

A second lens 1512 is arranged so as to image surface 1506 onto a detector 
1514, such as a CCD detector. Preferably the second lens 1512 is arranged such that the 
detector 1514 lies in its focal plane. The output of detector 1514, an example of which is 
a set of intensity maps designated by reference numeral 1515, is preferably supplied to 
data storage and processing circuitry 1516, which preferably carries out functionality 
"C" described hereinabove with reference to Fig. lA, providing an output indicating at 
least one and possibly both of the phase and the amplitude of the surface mapping 
wavefront. This output is preferably further processed to obtain information about the 
surface 1506, such as geometrical variations and reflectivity of the surface. 

In accordance with a preferred embodiment of the present invention, the beam 
of radiation supplied from radiation source 1500 has a narrow wavelength band about a 
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given central wavelength, causing the phase of the radiation reflected from surface 1506 
to be proportional to geometrical variations in the surface 1506, the proportion being an 
inverse linear function of the central wavelength of the radiation. 

In accordance with another preferred embodiment of the present invention, the 
beam of radiation supplied from radiation source 1500 has at least two narrow 
v/avelength bands, each centered about a different wavelength, designated Xu --j '^n- In 
such a case, the radiation reflected from the surface 1506 has at least two wavelength 
components, each centered around a wavelength A,i, and at least two indications 

of the phase of the surface mapping wavefront are obtained. Each such indication 
corresponds to a different wavelength component of the reflected radiation. These at 
least two indications may be subsequently combined to enable enhanced mapping of the 
surface 1506, by avoiding ambiguity in the mapping, knovm as 2n ambiguity, when the 
value of the mapping at a given spatial location in the surface exceeds the value of the 
mapping at a different spatial location in the surface by the largest of the different 
wavelengths X\, Xn^ A proper choice of the wavelengths X\, .... Xn , may lead to 
elimination of this ambiguity when the difference in values of the mapping at different 
locations is smaller than the muhiplication product of all the wavelengths. 

In accordance with still another preferred embodiment of the present invention, 
the phase manipulator 1510 applies a plurality of different spatial phase changes to the 
radiation wavefront reflected from surface 1506 and Fourier transformed by lens 1508. 
Application of the plurality of different spatial phase changes provides a plurality of 
differently phase changed transformed wavefronts which may be subsequently detected 
by detector 1514. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by phase manipulator 1510, 
resulting in at least three different intensity maps 1515. The at least three intensity maps 
are employed by the data storage and processing circuitry 1516 to obtain an output 
indicating at least the phase of the surface mapping wavefront. In such a case, the data 
storage and processing circuitry 1516, carries out functionality "C" described 
hereinabove with reference to Fig. lA, preferably in a maimer described hereinabove 
with reference to Fig. 13, where the wavefront being analyzed (Fig. 13) is the surface 
mapping wavefront. 
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Additionally, in accordance with a preferred embodiment of the present 
invention, the beam of radiation supplied from radiation source 1500 comprises a 
plurality of different wavelength components, thereby providing a plurality of 
wavelength components in the surface mapping wavefront and subsequently in the 
transformed wavefront impinging on phase manipulator 1510. In this case the phase 
manipulator may be an object, at least one of whose thickness, refractive index and 
surface geometry varies spatially. This spatial variance of the phase manipulator 
generates a different spatial phase change for each of the wavelength components, 
thereby providing a plurality of differently phase changed transformed wavefronts to be 
subsequently detected by detector 1514. 

Reference is now made to Fig. 16, which is a simplified partially schematic, 
partially pictorial illustration of a system for object inspection employing the 
functionahty and structure of Figs. lA and IB. As seen in Fig. 16, abeam of radiation, 
such as light or acoustic energy, is supplied from a radiation source 1600, optionally via 
a beam expander, onto at least partially transparent object to be inspected 1602. The 
radiation transmitted through the inspected object 1602, is an object inspection 
wavefront, which has an amplitude and a phase, and which contains information about 
the object 1602. At least part of the radiation transmitted through object 1602 is focused 
via a focusing lens 1604 onto a phase manipulator 1606, which is preferably located at 
the image plane of radiation source 1600. 

The phase manipulator 1606 may be, for example, a spatial light modulator or 
a series of different transparent, spatially non-uniform objects. It is appreciated that 
phase manipulator 1606 can be configured such that a substantial part of the radiation 
focused thereonto is reflected therefrom. Alternatively the phase manipulator 1606 can 
be configured such that a substantial part of the radiation focused thereonto is 
transmitted therethrough. 

A second lens 1608 is arranged so as to image object 1602 onto a detector 
1610, such as a CCD detector. Preferably, the second lens 1608 is arranged such that the 
detector 1610 lies in its focal plane. The output of detector 1610, an example of which is 
a set of intensity maps designated by reference numeral 1612, is preferably supplied to 
data storage and processing circuitry 1614, which preferably carries out functionality 
"C" described hereinabove with reference to Fig, lA, providing an output indicating at 



69 



least one and possibly both of the phase and the amplitude of the object inspection 
wavefront. This output is preferably further processed to obtain information about the 
object 1602, such as a mapping of the object's thickness, refractive index or 
transmission. 

In accordance with one preferred embodiment of the present invention, the 
beam of radiation supplied from radiation source 1600 has a narrow wavelength band 
about a given central wavelength, and the object 1602 is substantially uniform in 
material and other optical properties, causing the phase of the radiation transmitted 
through object 1602 to be proportional to thickness of the object 1602. 

In accordance with one more preferred embodiment of the present invention, 
the beam of radiation supplied from radiation source 1600 has a narrow wavelength 
band about a given central wavelength, and the object 1602 is substantially uniform in 
thickness, causing the phase of the radiation transmitted through object 1602 to be 
proportional to optical properties, such as refraction index or density, of the object 1602. 
It is appreciated that object 1602 may be any optical conduction element, such as an 
optical fiber. 

In accordance with another preferred embodiment of the present invention, the 
beam of radiation supplied from radiation source 1600 has at least two narrow 
wavelength bands, each centered about a different wavelength, designated In- In 

such a case, the radiation transmitted through object 1602 has at least two wavelength 
components, each centered around a wavelength Xi, and at least two indications 

of the phase of the object inspection wavefront are obtained. Each such indication 
corresponds to a different wavelength component of the transmitted radiation. These at 
least two indications may be subsequently combined to enable enhanced mapping of the 
properties, such as thickness, of object 1602, by avoiding ambiguity in the mapping, 
known as 2n ambiguity, when the value of the mapping at a given spatial location in the 
object exceeds the value of the mapping at a different spatial location in the object by 
the largest of the different wavelengths Xu - ^ ?^n* A proper choice of the wavelengths 
X,i, Xn , may lead to elimination of this ambiguity when the difference in values of 
the mapping at different locations is smaller than the multiplication product of all the 
wavelengths. 
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In accordance with still another preferred embodiment of the present invention, 
the phase manipulator 1606 applies a plurality of different spatial phase changes to the 
radiation wavefront transmitted through object 1602 and Fourier transformed by lens 
1604. Application of the plurality of different spatial phase changes produces a plurality 
of differently phase changed transformed wavefronts which may be subsequently 
detected by detector 1610. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by phase manipulator 1606, 
resulting in at least three different intensity maps 1612. The at least three intensity maps 
1612 are employed by the data storage and processing circuitry 1614 to obtain an output 
indicating at least the phase of the object inspection wavefront. In such a case, the data 
storage and processing circuitry 1614, carries out functionality "C" described 
hereinabove with reference to Fig. lA, preferably in a manner described hereinabove 
with reference to Fig. 13, where the wavefront being analyzed (Fig. 13) is the object 
inspection wavefront. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the beam of radiation supplied from radiation source 1600 comprises a 
plurality of different wavelength components, thereby providing a plurality of 
wavelength components in the object inspection wavefront and subsequently in the 
transformed wavefront impinging on phase manipulator 1606. In this case the phase 
manipulator 1606 may be an object, at least one of whose thickness, refractive index 
and surface geometry varies spatially. This spatial variance of the phase manipulator 
generates a different spatial phase change for each of the wavelength components, 
thereby providing a plurality of differently phase changed transformed wavefronts to be 
subsequently detected by detector 1610. 

Reference is now made to Fig. 17, which is a simplified partially schematic, 
partially pictorial illustration of a system for spectral analysis employing the 
functionality and structure of Figs. lA and IB. As seen in Fig. 17, a beam of radiation, 
such as light or acoustic energy, is supplied from a radiation source to be tested 1700, 
optionally via a beam expander, onto a known element 1702, such as an Etalon or a 
plurality of Etalons, Element 1702 is intended to generate an input wavefront, having at 
least varying phase or intensity. The radiation transmitted through the element 1702, is a 
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spectral analysis wavefront, which has an amplitude and a phase, and which contains 
information about the spectrum of the radiation source 1700, At least part of the 
radiation transmitted through element 1702 is focused via a focusing lens 1704 onto a 
phase manipulator 1706, which is preferably located at the image plane of radiation 
source 1700. 

The phase manipulator 1706 may be, for example, a spatial light modulator or 
a series of different transparent, spatially non-uniform objects. It is appreciated that 
phase manipulator 1706 can be configured such that a substantial part of the radiation 
focused thereonto is reflected therefrom. Alternatively the phase manipulator 1706 can 
be configured such that a substantial part of the radiation focused thereonto is 
transmitted therethrough. 

A second lens 1708 is arranged so as to image element 1702 onto a detector 
1710, such as a CCD detector. Preferably, the second lens 1708 is arranged such that the 
detector 1710 lies in its focal plane. The output of detector 1710, an example of which is 
a set of intensity maps designated by reference numeral 1712, is preferably supplied to 
data storage and processing circuitry 1714, which preferably carries out functionality 
"C" described hereinabove with reference to Fig. lA, providing an output indicating at 
least one and possibly both of the phase and the amplitude of the spectral analysis 
wavefront. This output is preferably further processed to obtain information about the 
radiation source 1700, such as the spectrum of the radiation supplied from radiation 
source 1700. 

In accordance with a preferred embodiment of the present invention, the 
spectral analysis wavefront is obtained by reflecting the radiation supplied from 
radiation source 1700 from element 1702, 

In accordance with another preferred embodiment of the present invention, the 
spectral analysis wavefront is obtained by transmitting the radiation supplied from 
radiation source 1700 through element 1702. 

In accordance with one more preferred embodiment of the present invention, 
the beam of radiation supplied from radiation source 1700 has a narrow wavelength 
band about a central wavelength, causing the phase of the radiation impinged on the 
object 1702 to be inversely proportional to the central wavelength supplied from 
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radiation source 1700 and related to at least one of a surface characteristic and thickness 
of element 1702. 

In accordance with another preferred embodiment of the present invention, the 
pluraHty of intensity maps 1712 are employed by the data storage and processing 
circuitry 1714, to obtain an output indicating at least one and possibly both of the phase 
and amplitude of the spectral analysis wavefront by expressing the plurality of intensity 
maps as at least one mathematical function of phase and amplitude of the spectral 
analysis wavefront and of plurality of different phase changes applied by phase 
manipulator 1706, wherein at least one and possibly both of the phase and amplitude is 
unknown and a function generating the different phase changes is known. This at least 
one mathematical function is subsequently employed to obtain an output indicating at 
least the phase of the spectral analysis wavefront. 

In accordance with still another preferred embodiment of the present invention, 
the phase manipulator 1 706 applies a plurality of different spatial phase changes to the 
radiation wavefront transmitted through element 1702 and Fourier transformed by lens 
1704, Application of the plurality of different spatial phase changes produces a plurality 
of differently phase changed transformed wavefronts which may be subsequently 
detected by detector 1710. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by phase manipulator 1706, 
resulting in at least three different intensity maps 1712. The at least three intensity maps 
are employed by the data storage and processing circuitry 1714 to obtain an output 
indicating at least the phase of the spectral analysis wavefront. In such a case, the data 
storage and processing circuitry 1714, carries out functionality "C" described 
hereinabove with reference to Fig. lA, preferably in a manner described hereinabove 
with reference to Fig. 13, where the wavefront being analyzed (Fig. 13) is the spectral 
analysis wavefront. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the beam of radiation supplied from radiation soiitce 1700 comprises a 
plurality of different wavelength components, thereby providing a plurality of 
wavelength components in the spectral analysis wavefront and subsequently in the 
transformed wavefront impinging on phase manipulator 1706. In this case the phase 
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manipulator may be an object, at least one of whose thickness, refractive index and 
surface geometry varies spatially. This spatial variance of the phase manipulator 
generates a different spatial phase change for each of the wavelength components, 
thereby providing a plurality of differently phase changed transformed wavefronts to be 
subsequently detected by detector 1710. 

In accordance with an embodiment of the present invention, the phase 
manipulator 1706 comprises a plurality of objects, each characterized in that at least one 
of its thickness and refractive index varies spatially. The spatial variance of the 
thickness or of the refractive index of the plurality of objects may be selected in a way 
such that the phase changes applied by phase manipulator 1706 differ to a selected 
predetermined extent for at least some of the wavelength components supplied by 
radiation source 1700. 

A specific selection of the objects is such that the phase change applied to an 
expected wavelength of radiation source differs substantially from the phase change 
applied to an actual wavelength of the radiation source. Alternatively, the spatial 
variance of the thickness or refractive index of the plurality of objects may be selected 
in a way such that the phase changes applied by phase manipulator 1706 are identical 
for at least some of the plurality of different wavelength components wavelength 
components supplied by radiation source 1700. 

In accordance with another embodiment of the present invention, the known 
element 1702 comprises a plurality of objects, each characterized in that at least one of 
its thickness and refractive index varies spatially. The spatial variance of the thickness 
or of the refractive index of the plurality of objects may be selected in a way such that 
the wavelength components of the input wavefront, generated by passing the 
wavelength components of the radiation supplied by radiation source 1700 through the 
element 1702, differ to a selected predetermined extent for at least some of the 
wavelength components supplied by radiation source 1700. 

A specific selection of the objects is such that the wavelength component of the 
input wavefront generated by an expected wavelength of radiation source differs 
substantially from the wavelength component of the input wavefront generated by an 
actual wavelength of the radiation source. Alternatively, the spatial variance of the 
thickness or refractive index of the plurality of objects may be selected in a way such 
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that the wavelength components of the input wavefront, generated by passing the 
wavelength components of the radiation supplied by radiation source 1700 through the 
element 1702, are identical for at least some of the wavelength components supplied by 
radiation source 1 700. 

Reference is now made to Fig. 18, which is a simplified partially schematic, 
partially pictorial illustration of a system for phase-change analysis employing the 
functionality and structure of Figs. lA and IB. As seen in Fig. 18, a known wavefront 
1800, which is a phase change analysis wavefront, having an amplitude and a phase, is 
focused via a focusing lens 1802, preferably performing a Fourier transform to 
wavefront 1800, onto a phase manipulator 1804, which is preferably located at the focal 
plane of lens 1802. The phase manipulator applies a plurality of different phase changes 
to the transformed phase change analysis wavefront. 

The phase manipulator 1804 may be, for example, a spatial light modulator or 
a series of different transparent, spatially non-uniform objects. It is appreciated that 
phase manipulator 1804 can be configured such that a substantial part of the radiation 
focused thereonto is reflected therefrom. Alternatively the phase manipulator 1804 can 
be configured such that a substantial part of the radiation focused thereonto is 
transmitted therethrough. 

A second lens 1806 is arranged so as to image wavefront 1800 onto a detector 
1808, such as a CCD detector. Preferably, the second lens 1806 is arranged such that the 
detector 1 808 lies in its focal plane. The output of detector 1808, an example of which is 
a set of intensity maps designated by reference numeral 1810, is preferably supplied to 
data storage and processing circuitry 1812, which employs the plurality of intensity 
maps to obtain an output indication of differences between the plurality of different 
phase changes applied by the phase manipulator 1804. 

In accordance with one preferred embodiment of the present invention, lateral 
shifts appear in the plurality of different phase changes. These may be produced, for 
example, by vibrations of the phase manipulator or by impurities in the phase 
manipulator. Consequently, corresponding changes appear in the plurality of intensity 
maps 1810, and result in obtaining an indication of these lateral shifts. 

In accordance with another preferred embodiment of the present invention, the 
plurality of intensity maps 1810 are employed by the data storage and processing 



75 



circuitry 1812 to obtain an output indicating the differences between the plurality of 
different phase changes applied by the phase manipulator 1804, by expressing the 
plurality of intensity maps as at least one mathematical function of phase and amplitude 
of the phase change analysis wavefront and of the plurality of different phase changes 
applied by phase manipulator 1804, where at least the phase and amplitude of the 
wavefront 1800 are known and the plurality of different phase changes are unknown. 
This at least one mathematical function is subsequently employed to obtain an output 
indicating at least the differences between the plurality of different phase changes. 

In accordance with still another preferred embodiment of the present invention, 
the phase manipulator 1804 applies a plurality of different spatial phase changes to the 
wavefront 1800 Fourier transformed by lens 1802, Application of the plurality of 
different spatial phase changes provides a plurality of differently phase changed 
transformed wavefronts which may be subsequently detected by detector 1808. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by phase manipulator 1804, 
resulting in at least three different intensity maps 1810. The at least three intensity maps 
are employed by the data storage and processing circuitry 1812 to obtain an output 
indicating at least the differences between the plurality of different phase changes. In 
such a case, the data storage and processing circuitry 1814, carries out functionality "C" 
described hereinabove with reference to Fig. lA, preferably in a manner similar to the 
manner described hereinabove with reference to Fig. 13, where the wavefront being 
analyzed (Fig. 13) is the known phase change analysis wavefront, and the spatial phase 
changes applied by phase manipulator 1804 are unknown. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the wavefront 1800 comprises a plurality of different wavelength 
components, thereby providing a plurality of wavelength components in the transformed 
wavefront impinging on phase manipulator 1804. In this case the phase manipulator 
may be an object, at least one of whose thickness, refractive index and surface geometry 
varies spatially. This spatial variance of the phase manipulator generates a different 
spatial phase change for each of the wavelength components, thereby providing a 
plurality of differently phase changed transformed wavefronts to be subsequently 
detected by detector 1 808. 
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Additionally, in accordance with another embodiment of the present invention, 
phase manipulator 1804 apphes one phase change to the radiation focused onto each 
spatial location thereon, resulting in one intensity map 1810, as an output of detector 
1808. In such a case, the data storage and processing circuitry 1812 employs the 
intensity map and the known wavefront 1800 to obtain at least an output indicating the 
phase change applied by phase manipulator 1804. 

In accordance with the foregoing methodology, the phase change applied by 
the phase manipulator may be a phase delay, having a value selected from one of a 
plurality of pre-determined values, including a possible value of zero phase delay and 
the output indication of the phase change obtained by data storage and processing 
circuitry 1812 is the value of the phase delay. In such a case, the phase manipulator may 
be media which stores information by different values of the phase delays at each of a 
multiplicity of different locations thereon, where the value of the phase delay constitutes 
the stored information. The stored information, encoded in the different values of the 
phase delays, is retrieved by data storage and processing circuitry 1812. It is appreciated 
that in such a case, wavefront 1800 may also comprise a plurality of different 
wavelength components, resulting in a plurality of intensity maps and consequently in 
an increase of the information encoded on the phase manipulator at each of the 
multiplicity of different locations. 

Reference is now made to Fig. 19, which is a simplified partially schematic, 
partially pictorial illustration of a system for stored data retrieval employing the 
functionality and structure of Figs. lA and IB. As seen in Fig. 19, optical storage media 
1900, such as a DVD or compact disk, has information encoded thereon by selecting the 
height of the media at each of a multipUcity of different locations thereon, as shown in 
enlargement and designated by reference numeral 1902. At each location on the media, 
the height of the media can be one of several given heights or levels. The specific level 
of the media at that location determines the information stored at that location. 

A beam of radiation, such as light or acoustic energy, is supplied from a 
radiation source 1904, such as a laser or a LED, optionally via a beam expander, onto a 
beam splitter 1906, which reflects at least part of the radiation onto the surface of the 
media 1900. The radiation reflected from an area 1908 on the media, onto which the 
radiation impinges, is a stored data retrieval wavefront, which has an amplitude and a 
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phase, and which contains information stored in area 1908. At least part of the radiation 
incident on area 1908 is reflected from the area 1908 and transmitted via the beam 
splitter 1906 onto an imaging system 1910, which may include a phase manipulator or 
other device which generates a varying phase function. 

Imaging system 1910 preferably carries out functionalities "A" and "B" 
described hereinabove with reference to Fig. lA, obtaining a plurality of differently 
phase changed transformed wavefronts corresponding to the stored data retrieval 
wavefront and obtaining a plurality of intensity maps of the plurality of phase changed 
transformed wavefronts. 

Preferably, imaging system 1910 comprises a first lens 1508 (Fig. 15), a phase 
manipulator 1510 (Fig. 15), a second lens 1512 (Fig. 15) and a detector 1514 (Fig. 15). 
The outputs of imaging system 1910 are supplied to data storage and processing 
circuitry 1912, which preferably carries out functionality "C" described hereinabove 
with reference to Fig. 1 A, providing an output indicating at least one and possibly both 
of the phase and amplitude of the stored data retrieval wavefront. This output is 
preferably further processed to read out the information encoded in area 1908 of media 
1900 and displayed on display unit 1914, 

In accordance with a preferred embodiment of the present invention, the beam 
of radiation supplied from radiation source 1904 has a narrow wavelength band about a 
given central wavelength, causing the phase of the radiation reflected from media 1 900 
to be proportional to geometrical variations in the media 1900, which contain the 
encoded information, the proportion being an inverse linear function of the central 
wavelength of the radiation. 

In accordance with another preferred embodiment of the present invention, the 
beam of radiation supphed from radiation source 1904 has at least two narrow 
wavelength bands, each centered about a different wavelength, designated X\, X^, In 
such a case, the radiation reflected from area 1908 in media 1900 has at least two 
wavelength components, each centered aroimd a wavelength Xi, >.n. 

At least two indications of the phase of the stored data retrieval wavefront are 
obtained, each such indication corresponding to a different wavelength component of 
the reflected radiation. These at least two indications may be subsequently combined to 
enhance mapping of the surface of area 1908 of media 1900 and therefore enhance 
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retrieval of the information, by avoiding an ambiguity in the mapping, known as 2% 
ambiguity, when the value of the height of the media at a given location exceeds the 
largest of the different wavelengths X\, , . Xn* 

In such a case, the range of possible heights at each location in area 1908 can 
exceed the value of the largest of the different wavelengths, without ambiguity in the 
reading of the heights. This extended dynamic range enables storing more information 
on media 1900 than would otherwise be possible. A proper choice of the wavelengths 
2.1, , may lead to elimination of this ambiguity when the difference of heights of 

the media in area 1908 at different locations is smaller than the multiplication product of 
all the wavelengths. 

In accordance with still another preferred embodiment of the present invention, 
a phase manipulator incorporated in imaging system 1910 applies a plurality of different 
spatial phase changes to the radiation wavefront reflected from media 1 900 and Fourier 
transformed by a lens, also incorporated in imaging system 1910. Application of the 
plurality of different spatial phase changes provides a plurality of differently phase 
changed transformed wavefronts which may be subsequently detected by a detector 
incorporated in imaging system 1910. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by a phase manipulator 
incorporated in imaging system 1910, resulting in an output from imaging system 1910 
of at least three different intensity maps. The at least three intensity maps are employed 
by the data storage and processing circuitry 1912 to obtain an output indicating at least 
the phase of the stored data retrieval wavefront. In such a case, the data storage and 
processing circuitry 1912, carries out functionality "C" described hereinabove with 
reference to Fig. lA, preferably in a manner described hereinabove with reference to 
Fig. 13, where the wavefront being analyzed (Fig. 13) is the stored data retrieval 
wavefront. 

Additionally, in accordance with an embodiment of the present invention, the 
beam of radiation supplied from radiation source 1904 comprises a plurality of different 
wavelength components, thereby providing a plurality of wavelength components in the 
stored data retrieval wavefront and subsequently in the transformed wavefront 
impinging on a phase manipulator incorporated into imaging system 1910, In this case 
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the phase manipulator may be an object, at least one of whose thickness, refractive 
index and surface geometry varies spatially. This spatial variance of the phase 
manipulator generates a different spatial phase change for each of the wavelength 
components, thereby providing a plurality of differently phase changed transformed 
wavefronts to be subsequently detected by a detector incorporated in imaging system 
1910, 

In accordance with another embodiment of the present invention, information 
is encoded on media 1 900 by selecting the height of the media at each given location to 
be such that the intensity value of the intensity map resulting from light reflected from 
the location and passing through imaging system 1910 lies within a predetermined range 
of values. This range corresponds to an element of the information stored at the 
location. By employing the plurality of intensity maps, multiple intensity values are 
realized for each location, each intensity value lying within a specific range of values. 
The resulting plurality of ranges of intensity values provide multiple elements of 
information for each location on the media 1900, 

It is appreciated that in such a case, retrieving the information stored at area 
1908 on the media from the outputs of imaging system 1910 may be performed by data 
storage and processing circuitry 1912 in a straight-forward manner, as by mapping each 
of the resulting intensity values at each location to their corresponding ranges, and 
subsequently to the information stored at the location. 

Preferably, the foregoing methodology also includes use of a phase 
manipulator incorporated in imaging system 1910, that appHes an at least time- varying 
phase change function to the transformed data retrieval wavefront impinging thereon. 
This time-varying phase change function provides the plurality of intensity maps. 

Alternatively or additionally, the beam of radiation supplied from radiation 
source 1904 comprises a plurality of different wavelength components, thereby 
providing a plurality of wavelength components in the stored data retrieval wavefront. 
The plurality of differently phase changed transformed wavefronts are obtained in 
imaging system 1910 by applying at least one phase change to the plurality of different 
wavelength components of the stored data retrieval wavefront. The phase changed 
transformed stored data retrieval wavefront can be detected by a single detector or 
alternatively separated, as by a dispersion element, into its separate plurality of different 
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wavelength components, each component being detected by a different detector. 

In accordance with yet another embodiment of the present invention, media 
1900 may have different reflectivity coefficients for the radiation supplied from light 
source 1904 at each of a multiplicity of different locations on the media. At each 
location on the media, a different percentage of the radiation may be reflected. The 
percentage may have one of several given values, where the specific value may at least 
partially determine the information stored at that location. 

In such a case, the information encoded on media 1900 may be encoded by 
selecting the height of the media at each of a multiplicity of different locations on the 
media and by selecting the reflectivity of the media at each of a multiplicity of different 
locations on the media. In such a case, more information can be stored at each location 
on the media, than could otherwise be stored. Furthermore, in such a case, employing an 
indication of the amplitude and phase of the stored data retrieval wavefront to obtain the 
encoded information includes employing the indication of the phase to obtain the 
information encoded by selecting the height of the media and employing the indication 
of the amplitude to obtain said information encoded by selecting the reflectivity. 

In accordance with still another embodiment of the present invention, the 
information is encoded onto media 1900 at several layers in the media. The information 
is encoded on the media by selecting the height of the media at each of multiplicity of 
different locations on each layer of the media. Each of these layers, placed one on top of 
the other in media 1900, is partially reflecting and partially transmitting. The beam of 
radiation from source 1904 impinging onto media 1900 is partially reflected from the 
top, first layer of the media, and partially transmitted to the layers lying therebelow. The 
energy transmitted by the second layer is partially reflected and partially transmitted to 
the layers lying therebelow, and so on, until the radiation tramsmitted through all the 
layers is partially reflected from the undermost layer. In such a case, radiation source 
1904 preferably includes a focusing system that focuses the radiation onto each one of 
the layers of media 1900 in order to retrieve the information stored on that layer. 
Alternatively or additionally, imaging system 1910 may include confocal microscopy 
elements operative to differentiate between the different layers. 

It is appreciated that area 1908 of media 1900 may be a relatively small area, 
comprising a single location on which information is encoded and possibly also 
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neighboring locations. In such a case, the detector incorporated in imaging system 1910 
may define only a single or several detection pixels. Additionally, the output indicating 
at least one and possibly both of the phase and amplitude of the stored data retrieval 
wavefront provided by circuitry 1912, includes at least one and possibly both of the 
height of the media and the reflectivity of the media at the location or locations covered 
by area 1908. 

In accordance with yet another embodiment of the present invention, the stored 
data retrieval wavefront comprises at least one one-dimensional component, 
corresponding to at least one one-dimensional part of area 1908 on media 1900, In such 
a case, the imaging system 1910 is preferably similar to the imaging system described 
hereinabove with reference to Fig. lOB. It preferably includes a first lens, such as 
cylindrical lens 1086 (Fig. lOB). 

The first lens preferably produces a one-dimensional Fourier transform, 
performed along an axis extending perpendicularly to the direction of propagation of the 
data retrieval wavefront, thereby providing at least one one-dimensional component of 
the transformed wavefront in a dimension perpendicular to direction of propagation. 
The first lens, such as lens 1086, focuses the stored data retrieval wavefront onto a 
phase manipulator, such as a single axis displaceable phase manipulator 1087 (Fig. 
lOB), which is preferably located at the focal plane of lens 1086. The phase manipulator 
1087 applies a plurality of different phase changes to each of the at least one 
one-dimensional components of the transformed wavefront, thereby obtaining at least 
one one-dimensional component of the plurality of phase changed transformed 
wavefronts. 

Additionally the imaging system may include a second lens, such as cylindrical 
lens 1088 (Fig. lOB), arranged so as to image the at least one one-dimensional 
component of the stored data retrieval wavefront onto a detector 1089, such as a CCD 
detector. Additionally the plurality of intensity maps are employed by circuitry 1912 to 
obtain an output indicating al least one and possibly both of the amplitude and phase of 
the at least one one-dimensional component of the data retrieval wavefront. 

Additionally, in accordance with the foregoing methodology, and as described 
hereinabove with reference to Fig. lOB, the phase manipulator 1087 preferably 
comprises a multiple local phase delay element, such as a spatially non-uniform 
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transparent object, typically including several different phase delay regions, each 
arranged to apply a phase delay to one of the at least one one-dimensional component at 
a given position of the object along a phase manipulator axis, extending perpendicularly 
to the direction of propagation of the wavefront and perpendicular to the axis of the 
transform produced by lens 1086. 

In such a case, there is provided relative movement betv^een the imaging 
system 1910 and the media 1900 along the phase manipulator axis. This relative 
movement sequentially matches different phase delay regions with different wavefront 
components, corresponding to different parts of area 1908 on media 1900, such that 
preferably each wavefront component passes through each phase delay region of the 
phase manipulator. 

It is appreciated that the relative movement between the imaging system 1910 
and the at least one one-dimensional wavefront component can be obtained by the 
rotation of media 1900 about its axis, while the imaging system is not moving. 

It is a particular feature of this embodiment, that each of the at least one 
one-dimensional component of the wavefront is separately processed. Thus, each of the 
at least one one-dimensional wavefront component, corresponding to a one-dimensional 
part of area 1908, is focused by a separate portion of the first cylindrical lens of imaging 
system 1910, is imaged by a corresponding separate portion of the second cylindrical 
lens and passes through a distinct region of the phase manipulator. The images of each 
of the one-dimensional parts of area 1908 at the detector incorporated in imaging 
system 1910 are thus separate and distinct images. It is appreciated that these images 
may appear on separate detectors or on a monolithic detector. 

In accordance with an embodiment of the present invention, the transform 
applied to the stored data retrieval wavefront includes an additional Fourier transform. 
This additional Fourier transform may be performed by the first cylindrical lens of 
imaging system 1910 or by an additional lens and is operative to minimize cross-talk 
between different one-dimensional components of the wavefront. In such a case, 
preferably an additional transform, such as that provided by an additional lens adjacent 
to the second cylindrical lens, is applied to the phase changed transformed wavelront. In 
such a case, preferably a further transform is applied to the phase changed transformed 
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wavefront. This further transform may be performed by the second cylindrical lens of 
imaging system 1910 or by an additional lens. 

Reference is now made to Fig. 20, which is a simplified partially schematic, 
partially pictorial illustration of a system for 3 -dimensional imaging employing the 
functionality and structure of Figs. lA and IB. As seen in Fig. 20, a beam of radiation, 
such as light or acoustic energy, is supplied from a radiation source 2000, optionally via 
a beam expander, onto a beam splitter 2004, which reflects at least part of the radiation 
onto a 3-dimensional object 2006 to be imaged. The radiation reflected from the object 
2006, is a 3-dimensional imaging wavefront, which has an amplitude and a phase, and 
which contains information about the object 2006. At least part of the radiation incident 
on the surface of object 2006 is reflected from the object 2006 and transmitted via the 
beam splitter 2004 and focused via a focusing lens 2008 onto a phase manipulator 2010, 
which is preferably located at the image plane of radiation source 2000. 

The phase manipulator 2010 may be, for example, a spatial light modulator or 
a series of different transparent, spatially non-uniform objects. It is appreciated that 
phase manipulator 2010 can be configured such that a substantial part of the radiation 
focused thereonto is reflected therefrom. Alternatively the phase manipulator 2010 can 
be configured such that a substantial part of the radiation focused thereonto is 
transmitted therethrough. 

A second lens 2012 is arranged so as to image object 2006 onto a detector 
2014, such as a CCD detector. Preferably the second lens 2012 is arranged such that the 
detector 2014 lies in its focal plane. The output of detector 2014, an example of which is 
a set of intensity maps designated by reference numeral 2015, is preferably supplied to 
data storage and processing circuitry 2016, which preferably carries out functionality 
"C* described hereinabove with reference to Fig. 1 A, providing an output indicating at 
least one and possibly both of the phase and amplitude of the 3-dimensional imaging 
wavefront. This output is preferably further processed to obtain information about the 
object 2006, such as the 3-dimensional shape of the object 

In accordance with a preferred embodiment of the present invention, the beam 
of radiation supplied from radiation source 2000 has a narrow wavelength band about a 
given central wavelength, causing the phase of the radiation reflected from object 2006 
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to be proportional to geometrical variations in the surface 2006, the proportion being an 
inverse linear function of the central wavelength of the radiation. 

In accordance with another preferred embodiment of the present invention^ the 
beam of radiation supplied from radiation source 2000 has at least two narrow 
wavelength bands, each centered about a different wavelength, designated ^i, A-n- In 
such a case, the radiation reflected from the object 2006 has at least two wavelength 
components, each centered around a wavelength Xi, and at least two indications 

of the phase of the 3 -dimensional imaging wavefiront are obtained. Each such indication 
corresponds to a different wavelength component of the reflected radiation. These at 
least two indications may be subsequently combined to enable enhanced imaging of the 
object 2006, by avoiding 2n ambiguity in the 3-dimensional imaging. 

In accordance with still another preferred embodiment of the present invention, 
the phase manipulator 2010 applies a plurality of different spatial phase changes to the 
radiation wavefront reflected from surface 2006 and Fourier transformed by lens 2008. 
Application of the plurality of different spatial phase changes provides a plurality of 
differently phase changed transformed wavefronts which may be subsequently detected 
by detector 2014. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial phase changes are applied by phase manipulator 2010, 
resulting in at least three different intensity maps 2015, The at least three intensity maps 
are employed by the data storage and processing circuitry 2016 to obtain an output 
indicating at least the phase of the 3-dimensional imaging wavefront. In such a case, the 
data storage and processing circuitry 2016, carries out functionality "C" described 
hereinabove with reference to Fig. lA, preferably in a manner described hereinabove 
with reference to Fig. 13, where the wavefront being analyzed (Fig. 13) is the 
3-dimensional imaging wavefront. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the beam of radiation supplied from radiation source 2000 comprises a 
plurality of different wavelength components, thereby providing a plurality of 
wavelength components in the 3-dimensional imaging wavefront and subsequently in 
the transformed wavefront impinging on phase manipulator 2010. In this case the phase 
manipulator 2010 may be an object, at least one of whose thickness, refractive index 
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and surface geometry varies spatially. This spatial variance of the ,phase manipulator 
generates a different spatial phase change for each of the wavelength components, 
thereby providing a plurality of differently phase changed transformed wavefronts to be 
subsequentiy detected by detector 2014. 

Reference is now made to Fig. 21 A, which is a simplified partially schematic, 
partially pictorial illustration of wavefront analysis functionality operative in 
accordance with another preferred embodiment of the present invention. The 
functionality of Fig. 21 A can be summarized as including the following 
sub-functionalities: 

A. obtaining a plurality of differently amplitude changed transformed wavefronts 
corresponding to a wavefront being analyzed, which has an amplitude and a phase; 

B. obtaining a plurality of intensity maps of the plurality of amplitude changed 
transformed wavefronts; and 

C. employing the plurality of intensity maps to obtain an output indicating at least 
one and possibly both of the phase and the amplitude of the wavefront being analyzed. 

As seen in Fig. 21 A, the first sub-functionality, designated "A" may be realized 
by the following functionalities: 

A wavefront, which may be represented by a plurality of point sources of 
light, is generally designated by reference numeral 2100. Wavefront 2100 has a phase 
characteristic which is typically spatially non-uniform, shown as a solid line and 
indicated generally by reference numeral 2102. Wavefi:ont 2100 also has an amplitude 
characteristic which is typically spatially non-uniform, shown as a dashed line and 
indicated generally by reference numeral 2103. Such a wavefront may be obtained in a 
conventional manner by receiving light from any suitable object, such as by reading an 
optical disk, for example, a DVD or compact disk 2104. 

A principal purpose of the present invention is to measure the phase 
characteristic, such as that indicated by reference numeral 2102, which is not readily 
measured. Another purpose of the present invention is to measure the amplitude 
characteristic, such as that indicated by reference numeral 2103 in an enhanced manner. 
A further purpose of the present invention is to measure both the phase characteristic 
2102 and the amplitude characteristic 2103. While there exist various techniques for 
carrying out such measurements, the present invention provides a methodology which is 



86 



believed to be superior to those presently known, inter alia due to its relative 
insensitivity to noise. 

A transform, indicated here symbolically by reference numeral 2106, is 
applied to the wavefront being analyzed 2100, thereby to obtain a transformed 
wavefront. A preferred transform is a Fourier transform. The resulting transformed 
v^ave front is symbolically indicated by reference numeral 2108. 

A plurality of different amplitude changes, preferably spatial amplitude 
changes, represented by optical attenuation components 2110, 2112 and 2114 are 
applied to the transformed wavefront 2108, thereby to obtain a plurality of differently 
amplitude changed transformed wavefronts, represented by reference numerals 2120, 
2122 and 2124 respectively. It is appreciated that the illustrated difference between the 
individual ones of the plurality of differently amplitude changed transformed 
wavefronts is that portions of the transformed wavefront are attenuated differently 
relative to the remainder thereof. 

As seen in Fig. 21 A, the second sub-functionality, designated "B", may be 
realized by applying a transform, preferably a Fourier transform, to the plurality of 
differently amplitude changed transformed wavefronts. Alternatively, the 
sub-functionality B may be realized without the use of a Fourier transform, such as by 
propagation of the differently amplitude changed transformed wavefronts over an 
extended space. Finally, functionality B requires detection of the intensity 
characteristics of plurality of differently amplitude changed transformed wavefronts. 
The outputs of such detection are the intensity maps, examples of which are designated 
by reference numerals 2130, 2132 and 2134. 

As seen in Fig. 21 A, the third sub-functionality, designated "C" may be 
realized by the following fimctionalities: 

expressing, such as by employing a computer 2136, the plurality of intensity 
maps, such as maps 2130, 2132 and 2134, as at least one mathematical function of 
phase and amplitude of the wavefront being analyzed and of the plurality of different 
amplitude changes, wherein at least one and possibly both of the phase and the 
amplitude are unknown and the plurality of different amplitude changes, typically 
represented by optical attenuation components 2110, 2112 and 2114 applied to the 
transformed wavefront 2108, are known; and 
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employing, such as by means of the computer 2136, the at least one 
mathematical function to obtain an indication of at least one and possibly both of the 
phase and the amplitude of the wavefront being analyzed, here represented by the phase 
function designated by reference numeral 2138 and the amplitude function designated 
by reference numeral 2139, which, as can be seen, respectively represent the phase 
characteristics 2102 and the amplitude characteristics 2103 of the wavefront 2100. In 
this example, wavefront 2100 may represent the information contained in the compact 
disk or DVD 2104. 

In accordance with an embodiment of the present invention, the plurality of 
intensity maps comprises at least four intensity maps. In such a case, employing the 
plurality of intensity maps to obtain an output indicating at least the phase of the 
wavefront being analyzed includes employing a plurality of combinations, each of at 
least three of the plurality of intensity maps, to provide a plurality of indications at least 
of the phase of the wavefront being analyzed. 

Preferably, the methodology also includes employing the plurality of 
indications of at least the phase of the wavefront being analyzed to provide an enhanced 
indication at least of the phase of the wavefront being analyzed. 

Also in accordance with an embodiment of the present invention, the plurality 
of intensity maps comprises at least four intensity maps. In such a case, employing the 
plurality of intensity maps to obtain an output indicating at least the amplitude of the 
wavefront being analyzed includes employing a plurality of combinations, each of at 
least three of the plurality of intensity maps^ to provide a plurality of indications at least 
of the amplitude of the wavefront being analyzed. 

Preferably, the methodology also includes employing the plurality of 
indications of at least the amplitude of the wavefront being analyzed to provide an 
enhanced indication at least of the amplitude of the wavefront being analyzed. 

It is appreciated that in this manner, enhanced indications of both phase and 
amplitude of the wavefront may be obtained. 

In accordance with a preferred embodiment of the present invention, at least 
some of the plurality of indications of the amplitude and phase are at least second order 
indications of the amplitude and phase of the wavefront being analyzed. 

In accordance with one preferred embodiment of the present invention, the 
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plurality of intensity maps are employed to provide an analytical output indicating the 
amplitude and phase. 

Preferably, the amplitude changed transformed wavefronts are obtained by 
interference of the wavefront being analyzed along a common optical path. 

In accordance with another preferred embodiment of the present invention, the 
plurality of intensity maps are employed to obtain an output indicating the phase of the 
wavefront being analyzed, which is substantially free from halo and shading off 
distortions, which are characteristic of many of the existing 'phase-contrast' methods. 

In accordance with still another embodiment of the present invention, the 
plurality of intensity maps may be employed to obtain an output indicating the phase of 
the wavefront being analyzed by combining the plurality of intensity maps into a second 
plurality of combined intensity maps, the second plurality being less than the first 
plurality, obtaining at least an output indicative of the phase of the wavefront being 
analyzed from each of the second plurality of combined intensity maps and combining 
the outputs to provide an enhanced indication of the phase of the wavefront being 
analyzed. 

In accordance with yet another embodiment of the present invention, the 
plurality of intensity maps may be employed to obtain an output indicating amplitude of 
the wavefront being analyzed by combining the plurality of intensity maps into a second 
plurality of combined intensity maps, the second plurality being less than the first 
plurality, obtaining at least an output indicative of the ampUtude of the wavefront being 
analyzed from each of the second plurality of combined intensity maps and combining 
the outputs to provide an enhanced indication of the amplitude of the wavefront being 
analyzed. 

Additionally in accordance with a preferred embodiment of the present 
invention, the foregoing methodology may be employed for obtaining a plurality of 
differently amplitude changed transformed wavefronts corresponding to a wavefront 
being analyzed, obtaining a plurality of intensity maps of the plurality of amplitude 
changed transformed wavefronts and employing the plurality of intensity maps to obtain 
an output of an at least second order indication of phase of the wavefront being 
analyzed. 

Additionally or alternatively in accordance with a preferred embodiment of the 
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present invention, the foregoing methodology may be employed for obtaining a plurality 
of differently amplitude changed transformed wavefronts corresponding to a wavefront 
being analyzed, obtaining a plurality of intensity maps of the plurality of amplitude 
changed transformed wavefronts and employing the plurality of intensity maps to obtain 
an output of an at least second order indication of amplitude of the wavefront being 
analyzed. 

In accordance with yet another embodiment of the present invention, the 
obtaining of the plurality of differently amplitude changed transformed wavefronts 
comprises applying a transform to the wavefront being analyzed, thereby to obtain a 
transformed wavefront, and then applying a plurality of different phase and amplitude 
changes to the transformed wavefront, where each of these changes can be a phase 
change, an amplitude change or a combined phase and amplitude change, thereby to 
obtain a plurality of differently phase and amplitude changed transformed wavefronts. 

In accordance with yet another embodiment of the present invention, a 
wavefront being analyzed comprises at least two wavelength components. In such a 
case, obtaining a plurality of intensity maps also includes dividing the amplitude 
changed transformed wavefronts according to the at least two wavelength components 
in order to obtain at least two wavelength components of the amplitude changed 
transformed wavefronts and in order to obtain at least two sets of intensity maps, each 
set corresponding to a different one of the at least two wavelength components of the 
amplitude changed transformed wavefronts. 

Subsequently, the plurality of intensity maps are employed to provide an output 
indicating the amplitude and phase of the wavefront being analyzed by obtaining an 
output indicative of the phase of the wavefront being analyzed from each of the at least 
two sets of intensity maps and combining the outputs to provide an enhanced indication 
of phase of the wavefront being analyzed. In the enhanced indication, there is no 2n 
ambiguity once the value of the phase exceeds 2%, which conventionally results when 
detecting phase of a single wavelength wavefront. 

It is appreciated that the wavefront being analyzed may be an acoustic radiation 
wavefront. 
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It is also appreciated that the wavefront being analyzed may be an 
electromagnetic radiation wavefront, of any suitable wavelength, such as visible light, 
infrared, ultra-violet and X-ray radiation. 

It is further appreciated that wavefront 2100 may be represented by a relatively 
small number of point sources and defined over a relatively small spatial region. In such 
a case, the detection of the intensity characteristics of the plurality of differently 
amplitude changed transformed wavefironts may be performed by a detector comprising 
only a single detection pixel or several detection pixels. Additionally, the output 
indicating at least one and possibly both of the phase and amplitude of the wavefront 
being analyzed may be provided by computer 2136 in a straight-forward manner. 

In accordance with an embodiment of the present invention, the plurality of 
different amplitude changes 2110, 21 12 and 21 14, preferably spatial amplitude changes, 
are effected by applying a time-varying spatial amplitude change to part of the 
transformed wavefront 2108, 

In accordance with a preferred embodiment of the present invention, the 
plurality of different amplitude changes 21 10, 21 12 and 21 14 are effected by applying a 
spatially uniform, time-varying spatial amplitude change to part of the transformed 
wavefront 2108. 

In accordance with an embodiment of the present invention, each of the 
wavefront 2100 and the transformed wavefiront 2108 comprises a plurality of different 
wavelength components. In such a case, the plurality of different spatial amplitude 
changes may be effected by applying an amplitude change to each of the plurality of 
different wavelength components of the transformed wavefront. It is appreciated that the 
amplitude changes may be spatially different or that the amplitude may be differently 
attenuated for each different wavelength component. 

In accordance with another embodiment of the present invention, each of the 
wavefront 2100 and the transformed wavefront 2108 comprises a plurality of different 
polarization components. In such a case, the plurality of different spatial amplitude 
changes may be effected by applying an amplitude change to each of the plurality of 
different polarization components of the transformed wavefront. It is appreciated that 
the amplitude changes may be spatially different or that the amplitude may be 
differently attenuated for each different polarization component. 
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In accordance with another embodiment of the present invention, the transform 
2106 applied to the wavefront 2100 is a Fourier transform, the plurality of different 
spatial ampUtude changes comprise at least three different amplitude changes, effected 
by applying a spatially uniform, time-varying spatial amplitude attenuation to part of the 
transformed wavefront 2108, and the plurality of intensity maps 2130, 2132 and 2134 
comprises at least three intensity maps. In such a case, employing the plurality of 
intensity maps to obtain an output indicating the amplitude and phase of the wavefront 
being analyzed preferably includes: 

expressing the wavefront being analyzed 2100 as a first complex function 
which has an amplitude and phase identical to the amplitude and phase of the wavefront 
being analyzed; 

expressing the plurality of intensity maps as a function of the first 
complex function and of a spatial fiinction governing the spatially uniform, 
time-varying spatial amplitude change; 

defining a second complex function, having an absolute value and a 
phase, as a convolution of the first complex function and of a Fourier transform of the 
spatial function governing the spatially uniform, time-varying spatial amplitude 
attenuation; 

expressing each of the plurality of intensity maps as a third fimction of: 

the amplitude of the wavefiront being analyzed; 

the absolute value of the second complex function; 

a difference between the phase of the wavefront being analyzed and 
the phase of the second complex function; and 

a known amplitude attenuation produced by one of the at least three 
different amplitude changes, to each of which one of the at least three intensity maps 
corresponds; 

solving the third function to obtain the amplitude of the wavefront being 
analyzed, the absolute value of the second complex function and the difference between 
the phase of the wavefront being analyzed and the phase of the second complex 
function; 

solving the second complex function to obtain the phase of the second 
complex function; and 
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obtaining the phase of the wavefront being analyzed by adding the phase 
of the second complex function to the difference between the phase of the wavefront 
being analyzed and the phase of the second complex function. 

Reference is now made to Fig. 2 IB, which is a simplified partially schematic, 
partially block diagram illustration of a wavefront analysis system suitable for carrying 
out the functionality of Fig. 21 A in accordance with a preferred embodiment of the 
present invention. As seen in Fig. 2 IB, a wavefront, here designated by reference 
numeral 2150 is focused, as by a lens 2152, onto an amplitude attenuator 2154, which is 
preferably located at the focal plane of lens 2152. The amplitude attenuator 2154 
generates amplitude changes, such as amplitude attenuation, and may be, for example, a 
spatial light modulator or a series of different partially transparent objects. 

A second lens 2156 is arranged so as to image wavefront 2150 onto a detector 
2158, such as a CCD detector. Preferably the second lens 2156 is arranged such that the 
detector 2158 lies in its focal plane. The output of detector 2158 is preferably supplied 
to data storage and processing circuitry 2160, which preferably carries out functionality 
"C" described hereinabove with reference to Fig. 21A. 

Reference is now made to Fig. 22, which is a simplified partially schematic, 
partially pictorial illustration of a system for surface mapping employing the 
functionality and structure of Figs. 21 A and 2 IB. As seen in Fig. 22, a beam of 
radiation, such as light or acoustic energy, is supplied from a radiation source 2200, 
optionally via a beam expander 2202, onto a beam splitter 2204, which reflects at least 
part of the radiation onto a surface 2206 to be inspected. The radiation reflected firom 
the inspected surface, is a surface mapping wavefront, which has an amplitude and a 
phase, and which contains information about the surface 2206. At least part of the 
radiation incident on surface 2206 is reflected from the surface 2206 and transmitted via 
the beam splitter 2204 and focused via a focusing lens 2208 onto an amplitude 
attenuator 2210, which is preferably located at the image plane of radiation source 2200. 

The amplitude attenuator 2210 may be, for example, a spatial light modulator 
or a series of different partially transparent non-spatially uniform objects. It is 
appreciated that amplitude attenuator 2210 can be configured such that a substantial part 
of the radiation focused thereonto is reflected therefrom. Alternatively the amplitude 
attenuator 2210 can be configured such that a substantial part of the radiation focused 
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thereonto is transmitted therethrough. 

A second lens 2212 is arranged so as to image surface 2206 onto a detector 
2214, such as a CCD detector. Preferably the second lens 2212 is arranged such that the 
detector 2214 lies in its focal plane. The output of detector 2214, an example of which is 
a set of intensity maps designated by reference numeral 2215, is preferably supplied to 
data storage and processing circuitry 2216, which preferably carries out functionality 
"C" described hereinabove with reference to Fig. 21A, providing an output indicating at 
least one and possibly both of the phase and the amplitude of the surface mapping 
wavefront. This output is preferably further processed to obtain information about the 
surface 2206, such as geometrical variations and reflectivity of the surface. 

In accordance with a preferred embodiment of the present invention, the beam 
of radiation supplied from radiation source 2200 has a narrow wavelength band about a 
given central wavelength, causing the phase of the radiation reflected from surface 2206 
to be proportional to geometrical variations in the surface 2206, the proportion being an 
inverse linear function of the central wavelength of the radiation. 

In accordance with an embodiment of the present invention, the beam of 
radiation supplied from radiation source 2200 has at least two narrow wavelength 
bands, each centered about a different wavelength, designated ,,,,XnAn such a case, 
the radiation reflected from the surface 2206 has at least two wavelength components, 
each centered around a wavelength Xi, ...^Xn . 

At least two indications of the phase of the surface mapping wavefront are 
obtained. Each such indication corresponds to a different wavelength component of the 
reflected radiation. These at least two indications may be subsequently combined to 
enable enhanced mapping of the surface 2206, by avoiding ambiguity in the mapping, 
known as In ambiguity, when the value of the mapping at a given spatial location in the 
surface exceeds the value of the mapping at a different spatial location in the surface by 
the largest of the different wavelengths Xu ,,.,Xn^ A proper choice of the wavelengths 
Xi, In , may lead to elimination of this ambiguity when the difference in values of 
the mapping at different locations is smaller than the multiplication product of all the 
wavelengths. 

In accordance with a preferred embodiment of the present invention, the 
amplitude attenuator 2210 applies a plurality of different spatial amplitude changes to 
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the radiation wavefront reflected from surface 2206 and Fourier transformed by lens 
2208. Application of the plurality of different spatial amplitude changes provides a 
plurality of differently amplitude changed transformed wavej&onts which may be 
subsequently detected by detector 2214. 

In accordance with yet another preferred embodiment of the present invention, 
at least three different spatial amplitude changes are applied by amplitude attenuator 
2210, resulting in at least three different intensity maps 2215. The at least three intensity 
maps are employed by the data storage and processing circuitry 2216 to obtain an output 
indicating at least one and possibly both of the phase and amplitude of the surface 
mapping wavefront. In such a case, the data storage and processing circuitry 2216, 
carries out functionality "C" described hereinabove with reference to Fig. 21 A, where 
the wavefront being analyzed (Fig. 21 A) is the surface mapping wavefront. 

Additionally, in accordance with a preferred embodiment of the present 
invention, the beam of radiation supplied from radiation source 2200 comprises a 
plurality of different wavelength components, thereby providing a plurality of 
wavelength components in the surface mapping wavefront and subsequently in the 
transformed wavefront impinging on amplitude attenuator 2210, In this case the 
amplitude attenuator may be an object, at least one of whose reflection and transmission 
varies spatially. This spatial variance of the amplitude attenuator generates a different 
spatial amplitude change for each of the wavelength components, thereby providing a 
plurality of differently amplitude changed transformed wavefronts to be subsequently 
detected by detector 2214. It is appreciated that the amplitude attenuation generated by 
attenuator 2210 may be different for each of the different wavelength components. 

In accordance with an embodiment of the present invention, the surface 2206 is 
a surface of media in which information is encoded by selecting the height of the media 
at each of a multiplicity of different locations on the media. In such a case, the 
indications of the amplitude and phase of the surface mapping wavefront provided by 
data storage and processing circuitry 2216 are employed to obtain the information 
encoded on the media. 

It is appreciated that other applications, such as those described hereinabove 
with respect to Figs. 16-20 may also be provided in accordance with the present 
invention wherein amplitude attenuation is performed instead of phase manipulation. It 
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is further appreciated that all of the applications described hereinabove with reference to 
Figs. 15-20 may also be provided in accordance with the present invention wherein 
both amplitude attenuation and phase manipulation are performed. 

It will be appreciated by persons skilled in the art that the present invention is 
not limited by what has been particularly shown and described hereinabove. Rather the 
present invention includes both combinations and subcombinations of features 
described hereinabove as well as modifications and variations of such features which 
would occur to a person of ordinary skill in the art upon reading the foregoing 
description and which are not in the prior art. 
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